NRL  Memorandum  Report  3077 


j 


CO 

pH 

o 

<c 


Localization  and  Dynamics  Determinations 
over  Spherical  Surfaces 

Albert  A.  Gerlach 

Advanced  Projects  Group 
A caustics  Di  vis  ion 


June  1975 


_ D D C 

;i;rr. 'v' 

^ AUG  8 1975 

1 Ikowtnstii 

D 


NAVAL  RESEARCH  LABORATORY 
Washington,  D.C. 


> 


Approved  for  public  releMie;  distribution  unlimited. 


3ECUHI  rv  CLASSlPiCATlON  or  T»,H  PAGE  'Whin  Out  Knl.i.d) 


REPORT  DOCUMENTATION  PAGE 


1 1 £r.i 

ID  INSTRUCTIONS 


READ  I 

BEFORE  COMPLETING  FORM 


• «tPOHT NUMBER 


NRL  Memorandum  Report  3077 


2 GOVT  ACCESSION  NO 


I.  RECIPIENT'S  CATALOG  NUMBER 


< title  (mao  iubuit*)-  ^ 

LOCALIZATION  AND  pYNAMICS  DETERMINATIONS 
OVER  SPHERICAL  SURFACE8  * _ 

"■  M w - V * 


NR 


S.  TYPE  OP  PEPOPT  • PCPIOO  COVCPCO 

interim  report  on  a continuing 
NRL  problem. 


t.  PEPPOPMINO  OPO.  PEPOPT  NUMBER 


I.  contpact  op  dp  an  t suSilRe 


7.  AutMORllJ 

Albert 


A^Gerlach  j (J^f  H |v  ! MS’  - 2^'/  7,1 


li.  PPOOPAM  ELEMENT.  PPOiECT.  task 
APE  A « WOPK  UNIT  NUMBERS 


» PEPPOPMING  ORGANIZATION  NAME  ARO 

Naval  Research  Laboratory 
Washington,  D.C.  20375 


AOOPESS 


- JNRL  Problem. S01-61 

I S ) v*  v*  ARPA  Ordetv2275  j 


II  CO1'’  BOLLING  OPPICE  NAME  AND  AOOPESS 

Defense  Advanced  Research  Projects  Agency 
Arlington,  VA  22209 


Ha 


2.  report  OATt 

JuniflP75  j 


T5  MONITORING  AGENCY  name  A AODRESSf/l  dllltrtn  from  Cortrolllni  Olllc.) 


. f , 

A V 4» 

/ 


J 


rr  number  op  rages 

93 


IS.  SECURITY  CLASS.  (»!  thlm  npMI| 

Unclassified 


ISa  DECL  ASSIPICATIOn/dOPNGRAOING 

schedule 


is  distribution  statement  foi  iMi  R.ponj 
Approved  for  public  release;  distribution  unlimited. 


' / < 

\ 


)//??/ , ■■  l;  7 / - L L I 

aginMia.  ■ 


17.  DISTRIBUTION  STATEMENT  (of  iho  abtfracf  ooftod  In  Block  20.  It  dIHoroot  from  Roport) 


P C 

; ■:  r 

■:  ■ r“r  ■ 3 -B?5 - 


i r, 


uLS^v.1..:.  .ij  U Ld 


D 


18  SUPPLEMENTARY  NOTES 


It  KEY  WORDS  fConrinut  on  f«v*ri<  iitf«  Nn*ct««ary  antf  Identify  by  block  numboe) 


\ 


Position  (location) 
Tracking 
Naviga.  m 

.Motional  dynamics 


Surveillance 
Spherical  geometry 
Analysis 


20  xXs TRACT  (Conilnum  on  iimp  .Id.  II  nccpatwr  and  Id.nttty  k>>  Slack  nuaikarj 

^Algorithms  providing  the  location  and  dynamics  (course  and  speed)  of  a moving  vehicle  we 
derived  for  application  to  navigation  and  localization  systems  on  a global  basis.  The  basic  mea- 
sures, used  for  computing  the  location  and  vehicle  dynamics,  are  obtained  from  signal*  propa- 
gated over  great-circle  paths  between  the  vehicle  in  question  and  two  (or  mow)  pain  of  wcemw 
sensors.  The  measured  signal  parameters  are  the  time  register  (or  time  delay)  and  the  time  scale- 
factor  (or  doppler)  shift,  between  pairs  of  receiving  sensors. 

^\^(Con  tinuet) 


DO  ijSn'tj  1473  EOlTION  OP  I NOV  •»  I*  OSIOLETC 


S/N  0 102-014’  SMI  I 


SECURITY  CLASSIFICATION  op  Twit  nn  RK  N*«  biAWRW 


) . 


/ / 


V. 

r 

-•  CLAitlflct  TIOW  or  tNH  PAOtrlHiwi  Pm  Enfitd) 


.20.  Abstract  (Continued) 

•^Algorithms  and  graphs  are  presented  to  further  provide  the  resolution  in  localization  and 
vehicle  dynamics  achievable  for  any  given  system  geometry.  These  resolution  (or  error)  data  are 
given  in  terms  of  the  resolution  in  tin*  basic  measures  of  time  register  and  time  scale-factor  shift. 
The  central  idea  of  the  report  is  to  provide,  in  convenient  and  usable  form,  algorithms  and  data 
which  yield  on  a global-system  basis  the  location  and  the  dynamics  of  a transiting  vehicle  and  the 
confidence  limits  for  these  parameters. 

Kxamples  are  provided  for  both  electromagnetic  signal  propagations  and  underwater  acoustic 
propagations  to  illustrate  the  feasibility  and  broad  applicability  of  the  resulting  data.  Optimum 
and  near  optimum  system  geometries  arc  derived  for  both  a worldwide  global  navigation  system 
and  for  a more  limited  surveillance  area.  The  results  indicate  that  a high  degree  of  resolution  in 
lx>th  localization  and  dynamics  is  achievable  over  extremely  long  ranges,  provided  reliable 
propagation  conditions  exist.  An  interesting  aspect  of  the  spherical  geometry  is  that  it  provides 
a natural  propagation  convergence  zone  on  the  opposite  side  of  the  globe  from  the  signal  source. 
This  phenomena  results  in  localization  and  dynamics  resolutions  at  extremely  long  range*  wh.cn 
are  equivalent  to  the  resolutions  achievable  in  the  near  sensor  field. 

It  is  conclud  ,U  that  full  * automatized  signal -processing  implementations  for  the  precise 
navigation  or  localization  ar,-'  "•citing  of  air,  sea,  and  land  vehicles  on  a global  scale  is  realizable 
in  the  near  future 
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LOCALIZATION  AND  DYNAMICS  DETERMINATIONS 
OVER  SPHERICAL  SURFACES 


I 

INTRODUCTION  AND  SUMMARY 

1.0  INTRODUCTION 

Navigation  or  localization  of  a transiting  vehicle  over  a wide  expanse 
of  the  earth's  surface  has  been  a problem  of  interest  to  man  since  the 
beginning  of  recorded  history.  Over  the  years,  new  techniques  for  this 
purpose  have  been  employed  as  the  relevant  technology  advanced.  At  the 
present  time,  technology  is  sufficiently  advanced  to  permit  localization 
on  a global  scale  to  well  within  one  square  mile,  using  energy  sensors  (or 
transmitters)  far  removed  from  the  point  under  surveillance.  The  major 
problems,  today,  consist  not  in  seeking  new  technology  but,  in  determining 
how  to  most  effectively  utilize  the  available  technology,  and  in  estab- 
lishing and  improving  the  reliability  of  anv  advanced  navigation  or 
localization  and  tracking  system. 

In  the  subject  paper  we  shall  address  tl\e  navigation  or  localization 
and  t-acking  problem  from  the  standpoint  of  exploiting  the  signal  properties 
of  long  range  propagations  which  travel  over,  essentially,  great  circle 
paths  between  an  unknown  location  on  the  surface  of  a sphere  (presumably  the 
earth)  and  several  known  locations  on  the  same  sphere.  The  approach  to  be 
taken,  in  the  subject  analyses,  is  to  postulate  a source  emanating  signal 
energy  from  the  unknown  point  in  question,  and  assume  that  the  signal  is 
being  received  and  processed  at  several  receiving  sensors  whose  locations 
are  precisely  known.  Under  these  circumstances,  the  basic  information 
about  the  source  is  remote  from  the  unknown  point  in  question.  When 
approached  from  this  standpoint,  the  problem  is  usually  termed  localization. 
The  problem  is,  however,  fully  reversible.  That  is,  the  same  analysis  is 
applicable  to  the  situation  where  a common  source  signal  is  transmitted  from 
some  several  precisely  known  locations  and  received  at  the  unknown  point  in 
question.  In  this  latter  case,  the  information  is  available  at  the  point  in 
question  to  determine  its  own  location.  Under  such  a situation,  the  self- 
localization  determination  is  generally  termed  as  navi' ation.  Thus,  from  a 
purely  theoretical  standpoint,  the  problem  of  localization  and  navigation  can 

Not?  Mdnuwnpt  submitted  M«*>  29,  197.'> 
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be  considered  reversible.  And  the  results  of  the  analyses  to  be  presented 
in  this  paper  will  be  applicnole  to  either  the  remote  localization  and 
tracking  of  an  unknown  located  signal  sw-ce,  or  to  the  navigation  (self- 
localization) of  a vehicle  transitinc.  .he  surface  of  the  sphere.  The 
reversible  feature  comes  about  as  a result  of  interchanging  the  role  of 
transmitting  source  and  receiving  sensors. 

2.0  BASIC  SYSTEM  CONCEPT 

The  basic  navigation  or  localization  concept,  as  studied  in  the  subject 
paper,  consists  of  exploiting  the  signal  measures  of  time  register,  t,  and 
time  scale-factor  shift,  i,  between  several  known  signal  locations.  Con- 
sider, if  vou  will,  a signal  emanating  from  a moving  source  and  propagating 
over  great  circle  paths  to  several  fixed  receiving  stat<cnt.  In  principle, 
the  signal  received  at  the  various  fixed  stations  will  be  irentical  except 
for  two  possible  factors.  The  one  factor  is;  the  received  signals  may  be 
out  of  time  register  (translated  in  time  with  respect  to  each  other)  depend- 
ing on  the  difference  in  great  circle  range  (between  the  source  and  the 
several  receiving  stations)  and  the  velocity  of  signal  propagation  in  the 
medium.  Thus,  if  one  can  consider  that  the  average  signal  propagation 
velocity  is  constant  and  known,  the  measure  of  the  time  register  'time 
shift)  between  the  signals  received  at  two  receiving  stations  will  provide 
a measure  of  the  difference  in  great  circle  range  between  the  two  receiving 
stations  and  the  source.  And,  as  will  be  shown  later,  the  acquisition  o' 
two  such  time  register  measures  (from  two  pairs  of  receiving  stations)  is 
sufficient  (in  a practical  sense)  to  uniquely  locate  the  source  on  the 
surface  of  the  sphere. 

The  second  factor  (or  difference  in  me  temporal  structure  of  the 
signals  received  at  me  several  receiving  stations)  is  a difference  in 
time  scale-factor.  This  difference  in  signal  structure  will  result  when 
the  motional  range-rate  of  the  source,  (along  the  great  circle  rays 
emanating  from  the  source  to  the  receiving  stations)  is  different,  i’he 
resulting  time  scale-factor  difference  is  manifested  as  a slight  time 
compression  (or  expansion)  ot  the  one  signal  relative  to  the  other.  This 
slight  shift  in  time  scale  iactoi  results  in  what  is  commc’ly  termed  a 
"doppler"  shift,  when  t he  signal  in  question  is  very  narrow  banded  in  Its 
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power  spe  trum.  Thus,  if  one  measures  the  time  scnle-factor  (or  doppler) 
shift  between  the  signals  received  at  two  receiving  stations,  he  will 
have  a meccure  of  the  difference  in  range-rate  of  the  source  along  the 
two  rays  emanating  from  the  source  to  the  two  receiving  stations.  And, 
as  will  be  shown  later,  the  acquisition  of  two  such  time  scale-factor 
shift  measures  (from  two  pairs  of  receiving  stations)  is  sufficient  to 
uniquely  determine  the  source  dynamics  (heading  and  speed). 

The  above  concept,  for  determining  localization  and  dynamics  informa- 
tion on  the  transiting  vehicle,  is  realizable  and  has  been  implemented  in 
a number  of  specific  applications  over  generally  relatively  short  ranges. 

A variety  of  techniques  are  available  to  the  system  designer  for  acnieving 
the  measures  of  time  register,  T,  and  time  scale-factor  shift,  5,  for  a 
given  application.  Depending  on  the  characteristics  of  the  transmitted 
signal  and  the  methods  of  achieving  the  t and  6 measures,  there  will  exist 
an  inherert  resolution  in  the  measure  of  these  critical  parameters.  That 
is,  in  a "real-life"  situation,  there  will  always  exist  some  fundamental 
limit  on  the  accuracy  with  which  the  time-register  measure,  t,  and  the 
time  scale-factor  shift  measure,  A,  can  be  relied  upon.  And,  this  accuracy 
or  resolution)  in  the  measures  of  t and  A is  generally  known  to  a reason- 
able degree  of  precision.  It  is  the  effect  of  this  imprecision  (in  the 
measure  of  ' and  A)  on  the  accuracy  of  the  computed  source  localization 
and  dynamics  which  is  of  great  interest  to  the  system  designer. 

In  the  subject  pajer,  we  intend  to  exploit  the  use  of  time  register 
measures  and  time  scale-factor  shift  measures,  between  several  receiving 
sensors,  to  determine  the  location  and  target  dynamics  (heading  and  speed) 
of  an  unknown  signal  source  located  on  the  surface  of  a sphere.  In  addi- 
tion, we  intend,  further,  to  determine  the  accuracy  (or  resolution)  of  the 
resulting  location,  heading,  and  speed  in  terms  of  the  (presumably)  known 
errors  (or  standard  deviations)  in  the  measures  of  t and  A.  Such  will  be 
the  object  of  the  sections  to  follow. 

3 • 0 DISCUSSION  AND  SUMMARY 

The  material  to  be  presented  is  divided  into  four  main  sections  to 
cover  the  four  logical  steps  of;  (1)  source  localization,  (2)  source 
dynamics  determination,  (3),  localization  resolution,  and  (4)  dynamics 
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resolution.  The  four  main  sections  are  reasonably  autonomous,  and  can 
be  read  in  any  order.  However,  as  one  might  expect,  the  sections  on 
resolution  do  build,  somewhat,  on  the  results  from  the  earlier  analyses. 

The  subject  material  is  presented  in  a concise  logical  manner,  in 
order  to  provide  the  reader  with  a thorough  understanding  of  the  under- 
lying physical  and/or  geometric  principles  involved.  Lengthy  mathematical 
proofs  and  derivations  are  avoided,  so  as  not  to  detract  from  the  principal 
theme  of  the  paper.  The  central  idea  of  the  paper  is  to  provide,  in  con- 
venient and  usable  form,  algorithms  and  data  which  yield;  (1)  the  location, 
(2)  the  dynamics,  and  (3)  the  confidence  limits  for  these  parameters,  for 
a transiting  vehicle  over  a global  surface.  The  assumed  system  inputs  are; 
(1)  the  precise  locations  of  several  receiving  sensors,  (2)  the  measures 
of  the  signal  physics  parameters  t and  ft,  and  (3)  the  resolution  of  these 
signal  physics  parameters  (standard  deviation  of  the  t and  ft  measures). 

From  the  analysis  and  data  provided,  one  will  also  be  able  to  synthesize 
a system  geometry  which  will  be  optimum  for  navigation,  or  localization 
and  tracking,  in  a specified  region  of  the  globe. 

In  the  first  two  sections  (Sections  II  and  III),  the  algorithms  are 
developed  for;  first,  the  source  location,  and  then  the  source  dynemics 
(heading  and  speed).  These  are  given  In  terms  of  the  known  locations  of 
the  fixed  sensor  positions,  and  the  measured  values  of  time  register,  T, 
or  time  scale-factor  shift,  ft.  A number  of  graphs  are  presented  which 
display  the  significant  geometric  properties  of  the  problem.  In  this  way, 
the  reader  may  obtain  an  insight  into  the  way  in  which  the  system  geometry 
influences  the  computed  source  location  and  dynamics  results. 

In  the  later  two  sections  (Sections  IV  and  V),  algorithms  and  data 
are  developed  which  provide  the  measures  of  confidence  that  can  be  placed 
in  the  earlier  computations.  These  are  given  in  terms  of  the  system 
geometry  and  the  resolution  <or  standard  deviation)  in  the  measurer,  of  t 
.•  d *.  A number  of  graphs  are  provided;  whereby,  one  m,-;''  quickly  ascertain 
the  confidence  in  the  computed  source  location  and  dynamics.  In  effect, 
one  may  readily  construct  a parallelogram  about  the  computed  source  location 
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within  which,  the  source  will  be  known  to  exist  with  a prescribed  level  of 
confidence.  Also,  the  heading  and  speed  of  the  source  will  be  known  within 
prescribed  standard  deviation  limits. 

Two  interesting  sidelights  resulting  from  the  subject  analysis  are 
worthy  of  special  mention  in  this  brief  summary.  The  first  of  these  con- 
cerns the  fact  that  the  spherical  geometry  results  in  a natural  propagation 
convergence  zone  on  the  diametric  opposite  side  of  the  globe  from  the  source 
location.  That  is,  rays  emanating  from  the  source  and  following  great 
circle  paths,  will  first  diverge  and  then  subsequently  converge  as  the 
signal  propagates  around  the  globe.  This  leads  to  some  very  interesting 
results  from  the  standpoint  of  accuracy  in  navigation  or  localization  and 
dynamics  determinations.  Under  normal  circumstances,  one  would  expect  that; 
as  a vehicle  is  located  more  and  more  remote  from  any  measurement  sensor 
configuration,  the  precision  with  which  that  vehicle  can  be  located  and 
its  dynamics  determined  will  degrade.  And  such  is  the  case,  up  to  the 
point  where  the  vehicle  is  one-quarter  of  the  distance  around  the  globe 
from  the  sensor  configuration.  However,  beyond  this  point,  the  navigation 
or  localisation  and  dynamics  precision  will  start  to  improve  again.  This 
improvement  will  continue  until  the  vehicle  reaches  the  extreme  fcr-field, 
on  the  opposite  side  of  the  globe  from  the  sensor  configuration.  In  this 
extreme  far-field,  the  localization  and  dynamics  resolution  will  be  the 
same  as  if  the  vehicle  were  located  in  the  image  near-field  of  the  sensor 
configuration.  This  convergence  phenomenon  has  some  very  interesting 
practical  implications.  It  implies  that  if  we  can  achieve  reliable  propa- 
gation over  the  surface  of  the  globe,  we  can  very  accurately  determine  the 
location,  heading,  and  speed  of  a vehicle  on  the  far  side  of  the  globe 
from  a simple  configuration  of  transmitters  (in  the  case  of  navigation) 
or  receiving  sensors  (in  the  case  of  localization  and  dynamics).  Typical 
resolutions  in  localization  and  dynamics  for  several  sensor  configurations 
are  given  in  the  examples  discussed  in  Sections  IV,  6,0  and  V,  4,0. 

The  second  interesting  sidelight,  resulting  from  the  subject  investi- 
gation, concerns  the  design  of  an  optimum  sensor  (or  transmitter)  configuration 
for  a global  navigation  system.  It  turns  out  that  the  localization  and 
dynamics  resolution  can  be  minimized  for  a geometric  configuration  of  trans- 
mitters located  on  three  orthogonal  globe  axes  over  the  surface  of  the  sphere. 
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Thus,  if  reliable  propagation  can  be  achieved  over  the  surface  of  the 
globe,  optimum  resolution  in  the  navigation  parameters  is  realizable.  Also, 
over  any  given  region  of  the  globe's  surface,  a more  confined  geometry  of 
sensors  can  be  designed  which  will  result  in  a resolution  of  the  navigation 
parameters  which  is  near  optimum.  Thus,  over  a limited  area  (such  as  an 
ocean  basin)  the  geometric  spacing  of  the  sensors  (or  transmitters)  need 
not  be  greater  than  the  diameter  of  the  area  under  surveillance.  These 
optimum  and  near  optimum  configurations  are  discussed  in  some  detail  in 
Sections  IV,  7.0  and  V,  5.0. 

4.0  CONCLUDING  REMARKS 

The  global  localization  and  dynamics  study,  presented  in  this  report, 
represents  an  extension  of  techniques  successfully  employed  over  relatively 
short  distances,  assuming  planar  geometry.  The  results  are  relevant  to  a 
variety  of  applications  in  the  general  areas  of  navigation  and  localization 
and  tracking  from  remote  sensor  configurations.  The  basic  measures  are 
readily  derived  from  the  signal  structure  of  energy  transmissions  which 
propagate  over,  essentially,  great  circle  paths  around  the  surface  of  the 
globe.  As  such,  the  results  are  applicable  to  both  underwater  aceustlc 
propagations,  and  to  the  lower  frequency  electromagnetic  emissions  which 
follow  (or  are  ducted  around)  the  surface  of  the  earth.  Current  technology 
is  amenable  today,  to  realize  both;  (1)  the  time  register  and  time  scale- 
factor  measurement  aspects,  and  (2)  the  automatic  computational  aspect:, 
of  a precise  global  navigation  or  localization  and  tracking  system.  There 
remains  to  be  investigate  in  some  detail,  the  deleterious  effects  of 
certain  long  range  propagation  anomalies  which  may  degrade  system  per- 
formance over  that  predicted  under  more  idealized  conditions.  It  can  be 
concluded  that  fully  automatized  signal  processing  implementations  for  the 
precise  navigation  or  localization  and  tracking  of  air,  sea,  and  land 
vehicles  on  a global  scale  is  realizable  in  the  near  future. 
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II 


SOURCE  LOCALIZATION 
OVER  SPHERICAL  SURFACES 

In  this  section,  we  shall  develop  the  algorithms  for  the  localization 
of  an  emitter  , source  whose  location  on  a spherical  surface  (presumably  the 
surface  of  the  earth)  is  unknown.  In  the  analyses  to  follow,  it  shall  be 
assumed  that  the  source  is  emitting  energy  which  propagates  in  great  circle 
paths  over  the  spherical  surface,  and  which  is  detectable  at  two  (or  more), 
pairs  of  receiving  sensors  whose  positions  on  the  surface  of  the  sphere  are 
precisely  known.  For  convenience,  Mu'  spherical  surface  shall  be  assumed 
to  be  the  earth,  and  all  location  coordinates  shall  be  given  in  angular 
units  ol  North  Latitude  and  West  Longitude. 

1 . 0 PROBLEM  FORMULATION  FOR  A RECEIVING  SENSOR  PAIR 

Consider  that  a signal  source  is  located  on  che  surface  of  a sphere 
at  location  and  that  a pair  of  receiving  sensors  are  detecting 

the  source  signal  at  known  locations  ,0, ) and  (Ov,0y).,  The  geometry 
of  the  problem  is  illustrated  in  Figure  1. 

In  the  subject  problem,  it  will  be  'ssumed  that  the  two  sensor  stations, 
x and  y,  can  measure  the  time  register  ioi  timi  difference  in  arrival) 
between  the  source  signals  which  arrive  at  the  two  receiving  stations.. 

Under  these  circumstances,  and  with  a-priori  knowledge  of  the  signal 
propagation  velocity,  the  range  difference  ."Ry,  is  readily  determined  as, 

' R*  , = Ry  - K„  = C "y  , ( 1 - 1 ), 

where , 

£Ky , is  the  great  circle  range  difference  between  the  source  to 

the  v receiving  station  and  the  source  to  the  x receiving 
station 

K,  * R,  , is  the  great  circle  range  from  the  source  to  die  x 
receiving  station 

= K^v  is  tl.e  great  circle  range  from  die  source  to  the  v 
receiving  st  at  ion 

is  the  time  difference  in  signal  propagation  between  the  Rv 
path  and  the  K,  path 

c is  the  mean  signal  propagation  velocity  over  the  great 

circle  pat h s 


t 


FIG.  1 GEOMETRY  OF  TWO  SENSOR  PROBLEM 


Since,  now,  Rx  and  Ry  will  be  a function  of  the  two  unknown  variables, 
G„  and  t,  , a unique  solution  to  the  source  location  will  not  be  possible  by 
a simple  measure  of  ilR^'or  Ty,).  What  equation  (1-1''  provides  is  a locus 
1 of  points  (or  curve)  over  the  surface  of  the  sphere;  such  that,  a signal 

source  located  anywhere  on  the  curve  would  provide  the  measure  T,,,  That 
is,  the  indicated  curve  maps  out  the  locus  of  points  on  the  surface  of  the 
sphere  where  the  great  circle  range  difference,  'iRy , ir  constant,  (This 
curve  is  depicted  in  Fig.  1 as  a dashed  line.)  To  obtain  a unique  solution 
for  tin  source  locution,  two  such  pairs  of  receiving  sensor  stations  will 
be  required.  Each  pair  will  provide  a locus  of  points  (curve)  over  the 
surface  of  the  sphere,  indicative  of  the  time  register  measure.  The  loca- 
tion of  the  sourct  will  be  the  point  where  the  two  curves  intersect.  (NOTE- 
Since,  in  fact,  the  two  curves  will  intersect  at  two  points  on  the  surface 
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ot  ' (<re,  ty  does  exist.  However,  practical  considerations 
Ck..  Tally  be  employed  to  rule  out  one  of  the  two  points  In  a real-life 
i i 3tion.) 


1 . 1 Fundamental  Equations : 

The  equation  which  describes  the  relevant  curve  for  a given 
measure  of  ARyl  (between  the  two  sensor  locations  x and  y)  may  take  a vari- 
ety of  forms.  We  shall  now  derive  several  of  these  forms  which  provide  a 
reasonably  clear  insight  Into  the  physical  geometry  of  the  problem. 


0 

It  can  readily  be  shown  that  the  chord  or  straight  line  distance 
between  two  arbitrary  points,  and  on  the  surface  of  a sphere, 

may  be  expressed  as, 


■ 2p  ^ l([l-8ln018in0a-cos01co89aco8(01-0;,)] 

-2 p J sin2lj(  0.j  - 0a)cos3l{(0i  -0a  )+8lna%(0l+0a)  slnalt(tf1  -flL,) 

I * 


- 2p  sin  (1-2) 

where, 

Rt,3  Is  the  chord  length  between  points  (9i,$i)  and  (0s,$a) 
on  the  surface  of  the  sphere 

p is  the  radius  of  Che  sphere 

is  the  angle  subtending  the  two  vectors  which  emanate 
from  the  center  of  the  sphere  to  the  two  points  in 
question 


and  where  all  0's  and  all  0's  are  expressed  in  angular  units  of  Latitude 
and  Longitude  lespectively. 


For  convenience,  we  shall  define  a new  variable,  as  follow*.  Let, 


Xm(6i  A:9a»^a>  * %[l-*ln01sineB-co80tco«0sco8(^1-tfc)] 


•? 

I 

't 

■i 

I 

i 

l 


it 


- slnal|(61-0a)co8Ol|(0, -#a)+cos8\i( 0^4-00 )slna\(tfi-^j>) 

- sln8[^v,(0l^,;0ot0o)  | (1-3) 
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X 


* 


It  will  also  prove  convenient  to  define  two  parameters,  u and  X as 
follows.  Let, 


A R>  , 


r « 


and, 


cTy„ 


Xyx  « sin(UyXsin_1^yXy')  ■ sin(ciyX /2p) 


(1-4) 


Since  | AR, » | will,  necessarily,  be  less  than  or  equal  to  RyX,  the 
value  of  ii  will  range  over  •!  s u < 1 , 


From  (1-2),  it  is  evident  that  the  great  circle  range  between  any  two 
points,  (9^0^  and  (0a,02),  on  the  surface  of  the  sphere  will  be, 

^i»a ( • ®a  »^a)l"2p  sin-1  (ft^s/2p)»2p  sin  ^ \\2  (1*5) 

Using  the  above  relations  in  (1-1),  then,  gives, 


or, 


Yy  (0,0;0y  ,«y  )-Yx  (e,0;0x  ,«5x)-UyxYXy  (0X  ,0X  ;9y  ,0y) 

sin-^X^O^Oy  ,0y)  ]-sin-lLx)(0,0;0x  ,0,)]  - uy*§£ 

■ Uy  X sin- X,  y(  0X  ,0X  ; 9y  ,0y  ) J 


(1-6) 


(1-7) 


1.2  Polynomial  Equations  in  X: 

Although  the  equations  given  in  (1-6)  and  (1-7)  display  the 
geometry  of  the  problem  in  simple  form,  they  are  far  too  cumbersome  for 
numerical  analysis.  By  taking  the  sine  of  both  sides  of  equation  (1-7)  and 
using  the  parameter  X defined  in  (1-4),  the  relation  may  be  transformed  to. 


- V X,(l-V)'  - >** 


d-8) 


And  by  a further  process  of  squaring  and  collecting  terms,  the  resulting 
equation  can  be  reduced  to  the  following  polynomial  equation  In  X,  and  \ » 

(X*  -X«  )®-2X*i  (X«+X|r  -2X«Xr  ) + Xj»  - 0 (l-»> 

# 
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or, 


(Xy  _Xy  « )°-2(l-2Xy  , )X»Xr  + (X«-Xy*)a  - X J , 


(1-10) 


The  above  equations  are  still  implicit  functions  of  the  two  variables 
0 and  0,  but  they  are  in  a form  suitable  for  the  solution  of  0 as  an 
explicit  function  of  the  variable  0. 

2.0  PROBLEM  SOLUTION  FOR  A RECEIVING  SENSOR  PAIR 

The  solution  to  equations  (1-9)  or  (1-10)  is  rather  involved  and 
tedious;  so  that,  the  details  will  not  be  provided  here.  We  shall  be  con- 
tent to  provide  only  the  final  results  and  demonstrate  the  use  of  these 
results  in  some  specific  examples  which  illustrate  the  general  properties 
of  the  solution. 

2 . 1  General  ‘.lgori  thins : 

For  the  sake  of  simplicity  the  pertinent  relations  are  given  in 
a nested  sequence  of  equations  which  start  with  the  final  format  of  6 
expressed  as  a function  of  0,  and  work  backwards  to  the  basic  coefficients 
which  are  simple  functions  of  the  sensor  location  parameters  (0»,0, ) and 

(0y  ,0 y). 


2.1.1  Format  of  Solution;  The  general  solution  of  the  subject 
problem  may  be  expressed  as, 


tan 


-l 


3*y(0;X,3,) 


i Xy  y (0  ; Xy  * y (0  > Xy  X ) 


(2-1) 


where  the  angles  S and  0 are  expressed  in  angular  units  of  North  Latitude 
and  West  Longitude  respectively,  and  where  Xy*  is  defined  in  (1-4).  In 
the  subject  application, 

Xyy  - Bina^(9jt  -8y  )cob3)((0*  *0y  )+cosa^(9x+6y  )sina^(0ii  *0y  ) (2-2) 

2.1.2  Functions  of  Longitude:  The  chree  basic  functious  of  the 

angle  0,  given  in  (2-1)  are. 


U 


3«y(0;\y*)  ■ f » ( ) (Xy * ) sin30  + 2f|C( ny)  (Xy x )sln0co80 


fe  ( x,y  ) (Xy  X )cOS  0 - f o (Xy  x ) 


•4y(0;Xy,)  - &,(W)  (Xy,)sin0  + g0(x,y)  (X?,)cos0 


(2-3) 

(2-4) 


^ x y ( 0 5 Xy  x ) ■ h,  ( *,y ) (Xy  x ) sin30  + 2h,  0 ( ) (Xy  x ) sin0CO30 


+ hc(x,y)  (Xyx)coss0  + h8(x,y)  (Xy«  ) 


(2-5) 


All  of  the  coefficients  of  the  trigonometric  functions  of  0,  in  the  above 
expressions,  are  polynomial  functions  of  the  parameter,  X3x  • 

2.1.3  Polynomial  Function  Coefficients:  The  sets  of  polynomial 

coefficients  given  in  (2-3),  (2-4),  and  (2-5)  are, 

^ « ( w ) (^y  * ) ” A*  ( *y ) + ( x,y ) ^y  * 


fc(x,y)(Xyx)  * Ae(x,y)  + Be(*y)Xyx 


3 2 

f 1C  ( *,y)  (Xyx  ) ■ A*o(x,y)  + D(*y)*y* 


fo(X3,)  = X,3,(l-X?,) 


2 2 
8«(w)(Xyx)  * E.(«,y)  +F,(w)Xyx 


8e  ( x,y ) ( Xy  x ) ■ Ef!  ( *y ) + Fc  ( *y  ) Xy  x 


(2-6) 


(2-7) 


*Mx,y)(Xyx)  * \y)  + H«(  j^y  ) Xy  X - B,  ( v ) xjx 

h«(*,y)(Xyx)  * Ge(x,y)  +He(x,y)Xyx-  B o(wr)Xjx 
h«o(x,y  ) (Xyx  ) ■ G«a(x,y)  + Hie(x,y)Xyx  ” D(^y)Xyx 
'’tfw)  (^?»)  ■ J(*y)  "^(*,y)^?*  ’ ^y* 


(2-8) 
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where,  the  coefficients  of  the  above  polynomials  In  Xs  are  exclusive 
functions  of  tho  receiving  acnaor  location*,  ) ami  (%,</>■,)• 


2.1.4  Fixed  Parameter  Coefficients:  Prior  to  giving  the 

functional  relations  fox'  the  fixed  parameter  coefficients,  the  following 
fundamental  parameters  will  be  defined  where  all  angles  0 and  0 are  given 
in  North  Latitude  and  West  Longitude  coordinates.  Let. 


Mfnr)" 

N<  *y)- 

p«  ( w ) 
pc(*y> 


1 -sln0,  sin0y 


lf(sin0,  -8in0y ) 


' %(cos0*  sin0* -cos0y  sin0j  ) 


I %(COS0«COB0,  -CO80y  COS0y  I 


(2-9) 


The  coefficients  of  the  polynomial  functions  given  In  (2-6),  (2-7), 
and  (2-8),  are  then, 


A.(W)  - p?(*y> 

(2-10) 

Ac  ( *y  ) ■ Pc  ( *.r ) 

(2-11) 

A,  e ( v ) * p • < ) Pe<JW) 

(2-12) 

B,(W)  ■ cos0,  cos0y  sin0,  8in0x 

(2-13) 

Bo<*r)  **  cos0*cos0yco80icos0y 

(2-14) 

D(*,y)  ■ %COS0*  COsOy  sin(0*+0y  ) 

(2-15) 

E.  ( «,y  ) ■ H{  H,y  ) P»  ( *y  ) 

(2-16) 

E«  ( %y ) ■ N(  ) P,  ( my  ) 

(2-l» 

IS 


F»(*,y)“  ^(cosO,  sinOy  sln0,  +sin0,  cosO,  sin0y  ) (2-18) 

Fe(»,y)”  \l(cOS0,  sinOy  CO80,+sin0y  CO80yCO80y  ) (2-19) 

G*  ( *y ) a %,y ) A,  ( x,y ) +N(  *,y ) (2Pi  ( ) P»(  ) -N(  5 B,  ( w ) ) (2-20) 

G*C  *.y)  “ M<  *r ) Ao(  n,y)  +N(  *,y ) (2Pe  (*y)  F«{  *y)  -N(  *y)  Be(*y)  1 (2-21) 


G*e  ( »y)  “M(  *y ) A«0  ( *y)  +N{  «,y)  [P»  ( «,y)  F«  (*,y ) +pe  ( *,y)  F»  ( *,y ) "N(  *y)  D(  w)  I (2-22) 
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N*  ( W ) ■ F»  ( nr ) *,y  ) B«  ( *y  ) ■ A,  ( ^y  ) 

(2-23) 

'*«(*.»)  ■ pc  ( «,y ) 4M(  ny  ) Be  ( ) - Afl(ny) 

(2-24) 

H.e(W)  - F*<*,y)  p0(ny)+M(».y)D(^,)  - Afo()%,) 

(2-25) 

J(*.y)  * N?v) 

(2-26) 

K(«,y)  ■ M(*y)  + J(*r) 

(2-27) 

L(  *y ) " 1 + M(  n,y ) 

(2-28) 

The  above,  fixed  parameter  coefficients,  are  dependent  only  on  the 
location  of  the  fixed  sensor  receiving  stations,  and  independent  of  both 
the  source  location  and  the  measurement  parameter  X . For  any  set  of 
fixed  sensor  locations,  x and  y,  the  above  parameters  may  be  computed  in 
advance  and  stored  in  memory.  It  is  further  easy  to  verify  that,  if 
Q(*,y)  represents  any  of  the  parameters  given  in  (2-10)  through  (2-28),  then, 

QOwr)  * <?(»,*>  (2-29) 

That  is,  the  fixed  paiameter  values  are  symmetric  with  respect  to  the  loca- 
tion subscripts  x and  y.  As  a consequence,  the  subscripts,  x and  y,  may  be 
transposed  without  affecting  the  results. 
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With  preknowledge  of  a aet  of  receiving  station  locations,  the  fixed 
parameter  coefficients  can  be  computed  and  tabulated  (or  stored  in  computer 
memory)  for  any  given  sensor  geometry. 


2. 2 Sign  Determination 

The  function  defined  by  (2-1)  describes  a family  of  curves  for 
0 as  a function  of  0 in  angulax*  coordinates  of  North  Latitude  and  West 
Longitude.  The  family  of  curves  is  generated  by  reason  of  the  measurement 
parameter  X . Thus,  when  X is  fixed  at  some  specified  value,  equation  (2-1) 
is  expected  to  generate  a single  closed  curve  over  the  surface  of  the  sphere. 
But  it  is  evident  that  some  ambiguity  exists  as  a result  of  the  "I  sign"  in 
the  denominator  of  equation  (2-1).  Unfortunately,  this  ambiguity  cannot  be 
resolved  without  placing  some  restrictions  on  the  location  of  the  two 
receiving  sensor  stations.  This  fact  arises  as  a result  of  the  possible  bl- 
valued  mature  of  the  independent  variable  8 in  (2-1).  Consequently,  some 
judgment  will  be  required  on  the  part  of  the  user  of  the  subject  algorithms 
if  they  are  to  be  applied  in  any  physical  application. 

The  basic  property  to  look  for  (in  a test  for  bl-valuedness)  is  to 
determine  whether  or  not  the  function  V»y(0;Xa)  can  become  negative  for 
some  value  of  0,  after  the  value  of  X is  specified.  If  so,  then  bl-valuedness 
of  0 is  implied  for  those  values  of  0 where  V’iy(0;Xa)  is  positive.  If  not, 
then  0 will  be  a single  valued  function  of  0 (for  the  specified  X)  and  only 
either  the  "+  sign"  or  the  "-  sign"  need  be  used  in  (2-1),  with  the  selection 
of  X carrying  its  own  "sign".  To  determine  the  appropriate  "sign"  to  use, 
in  this  latter  case,  the  following  rule  can  be  employed. 


for  0„  < 0y  use  the  "+  sign" 
for  0y  < 0„  use  the  "-  sign" 


(2-30) 


In  the  former  case  (where  bi-valuedness  for  0 is  indicated),  both  the 
"+  sign"  and  the  "-  sign"  will  need  to  be  employed  to  describe  the  appro- 
priate curve.  The  physical  geometry  of  the  sensor  station  locations, 

(8„ ,0, ) and  ($,,0^,  will  dictate  over  what  range  of  the  variable  0 the 
appropriate  curve  will  occur  for  a given  choice  of  the  parameter  X . 


In  general,  when  j 0„-0y  | «|  0, -9y  | , the  subject  family  of  curves 
will  be  single  valued  over  much  of  the  range  of  X . On  the  other  hand, 
when  | 0«-0y  1 «|  0,  -0y  | , the  subject  family  of  curves  will  be  bi-valued 
over  much  of  the  range  of  X. 


2.3  Divisions  of  Symmetry 

Due  to  the  physical  geometry  of  the  subject  problem,  the  family 
of  curves  (over  the  surface  of  the  sphere)  described  by  (2-1),  can  be 
expected  to  exhibit  symmetry  about  three  mutually  orthogonal  planes  which 
pass  through  the  center  of  the  sphere.  These  planes  will  Intersect  the 
surface  of  the  sphere  in  great  circles.  One  such  great  circle  curve  will 
pass  through  the  two  sensor  locations  (G,  ,0,)  and  (Sy,^).  A second  great 
circle  of  symmetry  (orthogonal  to  the  first)  will  bisect  the  great  circle 
arc  between  the  two  sensor  locations.  The  third  great  circle  of  symmetry 
will,  of  course,  be  unique  in  that  it  is  orthogonal  to  the  other  two. 

The  equations  for  the  great  circles  of  symmetry  are  of  interest  and 
may  be  readily  be  derived,  in  the  case  of  the  first  two. 


2*3.1  Algorithm  for  the  Great  Circle  Arc  Passing  Through  Two 
Fixed  Points,  (0,,0,)  and  (8y,0y):.  It  can  readily  be  shown  that  the  great 
circle  curve  which  passes  through  two  points,  (8«,0, ) and  (0y,0y)  is 
described  by  the  following  algorithm, 


0 


tan 


tan  9,  sin(0-0,  )+tan9,  sln(0v  -0)  I 
1 sin(0y-0"  ) i 


(2-31) 


where  all  9’s  and  all  0's  are  expressed  in  North  Latitude  and  West 
Longitude  coordinates  respectively.  (Note  that  for  0-0,, 0*6,  and  for 
0“0y  , fr*  % ) . 


2.3.2  Algorithm  for  the  Great  Circle  Arc  which  Bisects  and  is 
Orthogonal  tc  the  Great  Circle  Arc  Joining  Two  Fixed  Points,  (8»,0, ) and 
(9,  ,0y )z  To  determine  the  subject  curve,  we  need  merely  to  employ  (2-1), 
letting  the  parameter  Xy,*0.  When  this  is  done  we  obtain,  simply, 
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(2-32) 


tan" 


COsQg  COS(0-fl1)-CO89r,iOS(0»  -0) 
sin9y  -sin9x 


Vhe  above  great  circle  curve  will  be  orthogonal  to  the  great  circle 
curve  defined  by  (2-31)  and  will  intersect  this  curve  at  a point  midway 
between  the  two  fixed  points  ( 9»  ) and  (dj,0r). 

3.0  PROPER* IES  OF  THE  SOLUTION 

To  study  the  properties  and  characteristics  of  the  family  of  curves 
derived  for  the  sensor  pair  configuration,  it  will  be  convenient  to  trans- 
form the  coordinate  axes  of  the  sphere  in  a manner  which  will  provide  a 
simple  analytical  solution.  It  should  be  evident  that  the  relative  geometry 
of  the  family  of  curves  remains  invarient  with  a simple  set  of  rotations  of 
the  axes  of  the  sphere,  and  only  the  analytical  form  of  the  curves  change. 
(NOTE:  For  notatlonal  simplicity,  the  subscripts  will  be  dropped  from  the 

parameters  u and  X in  this  section.  It  should  be  kept  in  mind,  however, 
that  the  indicated  subscripts  are  a proper  part  of  these  parameters). 

3.1  Longitudinal  Alignment  of  the  Receiving  Sensors: 

We  now  propose  a simple  set  of  rotations  on  the  spherical  axes; 
such  that,  the  transformed  coordinate  points  of  the  receiving  sensor  stations 
will  be  aligned  along  the  aero  Longitudinal  axis  and  equidistant  about  the 
equator.  That  is,  let 


0r  ■ -e, 

0,  - 0,  - 0 


(3-1) 


Under  these  circumstances,  the  equation  for  the  family  of  curves  given 
by  (2-1)  becomes. 


6 - tan"1  ^ 


(l-la)-C0Sae,C088fl  I 
. s.  7 a ' 


tan 


-i 


7TP 
x 

yi-x3  'y 


f — 2 /.  co8rt0*  sin*#  ) 

1 rprs  a / 1 + . a y 


sin  8,  • X > 

sina9,  -XJ  J 


(3-2) 
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Now,  from  ('-4)  and  (2-2),  we  may  write 


Xxy  = sln20x 
X = sinu0x 


Using  (3-3)  in  (3-2)  gives, 


0 = tan 


•r./i  +£££ 

V sin 


s8„  sins0 


0X  -sin  u8x 


tanu  0, 


tan 


,i  f L cos  8X  sin  0 


.tanu  9X 


sin(l+u)  flx  sin(l-u)  6X 


(3-3) 


(3-4) 


3.1.1  Case  When  9„  » n/2;  In  the  particular  case  where 
0j  - rt/2,  equation  (3-4)  reduces  to,  simply, 

6 - utt/2  (3-5) 

In  other  words,  the  subject  family  of  curves  degenerate  to  simple  Latitude 
circles  over  the  surface  of  the  sphere.  This  is  not  surprising;  since,  for 
receiving  sensors  located  at  the  two  poles,  the  circles  of  constant  Latitude 
represent  the  locus  of  points  (on  the  surface  of  the  sphere)  whose  difference 
in  great  circle  path  length  to  the  two  poles  is  a constant. 

3.1.2  Symaetrlcal  Properties:-  An  examination  of  equation  (3-2) 
reveals  that 

0(-0;X)  « C ( 0 ; X ) 

0(TT-0,\)=  0,'0;\)  (3-6) 

0(0 ; - X)  = -O(0;X) 

Therefore,  the  family  of  curves  defined  in  (3-4)  is  symmetric  about  each  of 
three  mutually  orthogonal  planes  which  pass  through  the  center  of  the  sphere. 

Ln  this  particular  coordinate  system,  the  three  planes  intersect  the  surface 
of  the  sphere  in  the  three  orthogonal  great  circles  which  form  the  zero  and 
180  degree  Longitudinal  curves,  the  ± 90  degree  Longitudinal  curves  and  the 
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zero  Latitude  curve  (equator).  This  symmetrical  property  of  the  subject 
family  of  curves  had  been  forecast  in  Section  2,3  . 

3.1.3  Maxima  and  Minima  of  1 9 j From  (3-4)  it  may  readily  be 
observed  that  the  minima  and  maxima  of  the  absolute  value  of  9(0;\)  will 
occur  for  0 equal  to  zero  (and  380  degrees)  and  ± n/2  (or  ± 90  degrees) 
respectively.  The  minimum  absolute  value  of  6 will,  therefore,  be 

6.  = U0,  (3-7) 


and  the  maximum  absolute  value  of  0 becomes, 


0M  ■ tan*1 


iinuO, 


in?0,  -sin2u0„ 


. sinu  0, 

tan*  1 


V«in  ( 1+u ) 0,  sin ( 1 -u ) fC 


(3-8) 


A 

3.1.4  Case  When  1(  When  0,  is  sufficiently  small  so  that 

sinG,  and  tanO,  can  be  replaced  by  6,  without  significant  error,  equations 
(3-4)  and  (3-8)  can  be  approximated  as, 


and 


0*.  ~ tan 


sin* 


(3-9) 


0 = tan*1 


cos‘;  A»  + sinH0  tan0M  ^ 


(3-10) 


The  above  system  will  now  be  examined  in  the  near  field  (relatively 
close  proximity  to  t tie  receiving  sensors),  in  the  relatively  far  field 
(remote  from  the  receiving  sensor  stations),  and  in  the  extreme  far  field 
(the  image,  of  the  near  field  on  the  opposite  side  of  the  sphere). 


When  0 and  0 are  sufficiently  small;  such  that,  tan  " and  sin  0 can 
be  replaced  by  the  respective  angles  without  significant  error,  then 
equation  (1-10)  can  be  approximated  as, 
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1 


(3-11) 


f 02 

’ o-u2)e,2 

The  above  relation  is  the  equation  for  e family  of  hyperbolas  in  the 
9,0  plane.  This  is  just  what  one  would  expect,  since  the  solution  for  the 
basic  phenomena  in  planar  geometry  results  in  a hyperbolic  family  of  curves. 
A graphical  plot  of  (3-11)  for  the  near  field  solution  is  illustrated  in 
Figure  2. 

The  relatively  far  field  behavior  of  the  family  of  curves  can  be 
obtained  from  (3-10)  through  the  artifice  of  taking  the  limit  as  0X-*  0. 

Thu  s , 


Lim.  9 - tan'1 
9,-0 


r u 


T sin0  | 


» tan*1  ftan0Msin0 ) 


(3-12) 


A plot  of  the  resulting  family  of  curves  is  illustrated  in  Figure  3. 

An  Interesting  facet  of  the  spherical  geometry  is  that,  in  the 
extreme  far  field  (the  surface  area  on  the  opposite  side  of  the  sphere  from 
the  receiving  sensor  location),  there  exists  a convergence  zone  where  the 
ray  paths  converge  to  form  a propagation  field  identical  to  the  field  in 
the  near-sensor  region.  As  r.  consequence,  the  equations  and  curves  given 
for  the  near  field  are  equally  applicable  to  the  extreme  far  field  region. 

This  fact  has  important  practical  implications;  in  that,  it  tells  us  that 
the  resolution  for  localizing  a source  (or,  conversely,  the  accuracy  of 
navigation)  in  the  extreme  far  field  (on  the  opposite  side  of  the  sphere) 
is  just  as  good  as  in  the  near  field.  This  is  a consequence  of  the  propa- 
gation ray  convergence  in  the  extreme  far  field  zone. 

3.1.5  Case  When  is  large:-  T.'  study  the  behavior  of  the  curves 
for  the  case  when  0,  is  relatively  large,  we  shall  let  0,  be  equal  to 
'Vb  (30  degrees).  In  this  event,  equation  (3-4)  becomes, 

0 - tan*1  L/l  + — CEpll tan(un/6)}  (3-13) 

V l-4sin‘ (uv/6)  J 


i 
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A graphical  presentation  of  these  curves  is  illustrated  in  Figure  4. 
The  family  of  curves  is  plotted  over  the  entire  range  of  Longitude,  and  for 
both  positive  and  negative  values  of  the  measurement  parameter  i.  The 
illustration  vividly  demonstrates  the  three  dimensions  of  symmetry  over  the 
surface  of  the  sphere.  The  location  of  the  two  receiving  sensors  are 
depicted  as  small  circles  located  at  6 = ± 30 0 and  0 = 0.  In  a sense,  the 
family  of  curve',  take  on  the  appearance  of  flexible  Latitude  circles  which 
have  been  pinched  (or  compressed)  toward  the  equator  (9=0),  at  0=0  and 
0=180°,  due  to  the  presence  of  the  two  sensors.  The  closer  the  two  sensors 
are  spaced,  the  greater  will  be  the  compression  of  the  family  of  curves  at 
the  two  relevant  Longitudinal  positions. 


3.1,6  Slope  of  the  Curves:'  To  determine  the  slope  of  the  family 
of  curves  given  by  (3-4)  we  may  proceed  as  follows, 


Slope 


1 d§ 
cosO  d0 


cos?0, sinu 0, sin0cos0 


. r 


'\j  (cos'  u 6,  -cos"  6,  cos‘0)  (sin?0,  cos  u0,  - cos'  0,  sin'uG,  cos'  0) 


(3-14) 


When  0«  is  small  this  approximates, 

wsin0cos0  (1-9;) 

Slope  = -y=zz====z  z======z======== : . ======^z======z~. - — — (3-15) 

aJ  [ l-u‘  (cos "0+0!  ' sin'  0)3 | sin"  0+9, ' (cos~0-u‘ ) ] 


And,  in  the  limit,  as  0,  -*  U, 

acos0 

Slope 


Vl-uacO8?0 


(3-16) 


Thus,  for  0 *<  1,  the  slope  becomes  equal  to  u/«/  l-u"  , which  is  the 
slope  asymptote  for  the  plane  geometry  case  (see  Figure  2). 

3 . 2 Fquatorial  Alignment  of  the  Receiving  Sensors 

As  an  alternate  to  aligning  the  receiving  sensors  in  a longitudinal 
orientation,  we  could  have  transformed  the  coordinate  axes  ot  the  sphere  to 
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TON  OF  $ AS  A FUNCTION  OF  $ for!  5,  BQUAL  TO  30° 


a new  set  where  the  receiving  sensors  were  located  on  the  e^ator,  equi- 
distant about  the  zero  Longitudinal  axis.  That  is,  let', 


6y  = 8,  = 0 

0y  = -<t>X 


Under  these  circumstances,  the  equation  for  the  family  of  curves 
given  by  (2-1)  becomes, 

-\ 


r>  , . 


, . J(ain20%~\2)  (gini!0-XB)  , 

0 = ± tan'1  < ta  / J 


(3-18) 


Now,  from  (1-4)  and  (2-2),  we  may  write, 


Xxy  s sin?0, 
X * sinu0. 


(3-19) 


Using  ( 3- 1 9 ) in  (3-18)  gives 


fs/(  sin  E0: 
* -t  tan"1 1 


, -sin“u0, ) (sinB0-sin‘  u0, ) 


sinu0*  cosu0. 


«=  ± tan‘ 


sin*  0,  _ 

Tin7!?! 


\ fsln‘0 

/ ' cns'^&t 


tan'  u0. 


(3-20) 


The  resulting  family  of  curves  given  by  the  above  expression  is 
identical  to  the  family  of  curves  given  in  (3-4).  The  only  difference 
between  the  two  is  a change  in  the  coordinate  frame  of  reference.  As  a 
consequence,  all  of  the  properties  developed  in  Section  3,1  will  be 
applicable  to  the  above  curves.  It  is  important  to  realize  that  the 
geometric  properties  of  the  family  of  curves  are  independent  of  any  linear 
transformation  on  the  frame  of  reference  over  the  surface  of  the  sphere. 
However,  it  should  be  noted  that  an  appropriate  linear  transformation  in 
the  co  rdinate  frame  of  reference  can  greatly  simplify  the  resulting 
analytical  representation  of  the  family  of  curves. 
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The  coordinate  frames  of  reference  chosen  in  the  last  two  sections 
provide  the  simplest  forms  for  the  analysis  of  the  geometric  properties  of 
the  subject  family  of  curves.  In  fact,  the  coordinate  frame  of  reference 
chosen  in  Section  3.1  Is  superior  to  this  later  choice  since  the  ambiguity 
relative  to  the  "sign"  of  the  parameter  u was  removed.  This  was  evident 
in  (3-4)  where  the  "±  sign"  was  not  present.  (The  parameter  u can  carry 
its  own  "sign"  when  0,  - 0V  ) . In  (3-20),  on  the  other  hand,  it  was 
necessary  to  incorporate  the  "±  sign"  in  order  to  completely  map  the  curve 
for  any  fixed  value  of  the  parameter  u.  In  this  latter  case,  one  would 
select  the  curve  for  positive  values  of  <t>  when  uyx  is  positive,  and  the 
curve  for  negative  values  of  0 when  uy „ is  negative.  In  the  more  general 
case,  given  by  (2-1),  some  care  will  need  to  be  taken  (In  this  regard)  if 
ambiguity  in  the  location  of  the  proper  curve  i3  to  be  avoided. 

4 . 0 SOURCE  ;.(x:alix.ation  algorithms 

In  the  earlier  analyses  we  have  determined  that  a simple  measure  of 
range  (or  time)  difference  between  a source  and  two  known  receiving  sensor 
stations  Is  insufficient,  b>  itself,  in  determining  the  location  of  the 
source.  Such  a measure  will,  however,  provide  a closed  curve  over  the 
•urface  of  the  sphere  on  which  the  source  is  known  to  be  located.  Conse- 
" ently,  if  the  range  (or  time)  difference  between  two  (or  more)  pairs  of 
k h -i,  receiving  sensor  locations  can  be  established,  then  the  location  r r 

3C..'Vo  will  be  known  to  exist  at  the  intersection  of  the  two  (or  mor. 
'■su.iiiu'  closed  curves.  Since  two  closed  curves  will  Intersect  at  two 
def  . minl'-ric  locations  on  the  surface  of  the  sphere,  bivalued  ambiguity 

exists  its  principle.  However,  one  can  generally  rule  out  one  of  the 
on  the  basis  oi  purely  physical  considerations.  As  a consequence, 
thi  ’ '.lization  of  an  unknown  source  is,  for  all  practical  purposes, 
real  i.-'io'-e  in  physical  applications. 

4.1  S irm, I taneous  Equations 

t,‘  l-  invoked  in  the  localization  problem.  Is  the  simultaneous 
solution  of  “-■)  .Mjuations  of  the  form  given  in  (2-1)  for  each  of  two  sensor 
pairs.  Thus,  _l;e  ocatiori  (in  North  I.atitude  and  West  Longitude  coordinates) 
of  the  source  wl.l  do  the  set  (0o ,^ ) yhich  satisfies  the  following  relations, 
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G0  = tan"1  ^ 


-7,  y (00  » )-y  x ) 


i V *<y  y (0o  > * )'  "^ty  (0o  i 'ly x ) 


tan 


.1 


7tw(0o;^r) 


± \*  z-J  « (0o ; X«  i ')~'&i » (0o ; X,  t ) 


(4-1) 


where,  the  known  locations  of  the  receiving  sensors  are  ( 0, ,0, ) , ( Gy ,0y ) , 

( Gj ,0, ) and  (0»,0V).  (It  should  be  aoted  that  only  three  receiving  sensor 
locations  are  necessary,  since  two  pairs  can  readily  be  derived  from  only 
three  sensor  stations.  In  this  event,  z can  be  made  equal  to  y in  the 
above  relations.) 

4.2  General  Approaches  to  the  Solution 

Now,  there  are  a number  of  ways  for  solving  the  above  relations. 
One  simple  technique  is  to  maintain  a graphical  plot  of  the  family  of 
curves  resulting  from  the  two  sensor  pairs.  Upon  obtaining  the  measures 
for  and  u%1,  one  may  then  identify  and  follow  the  trace  for  these  two 
parameters  on  the  graphical  plots  until  the  appropriate  point  of  inter- 
section ia  obtained.  Although  the  graphical  technique  is  rather  crude  and 
generally  not  too  accurate  (unless  extremely  large  graphs  are  maintained 
with  quite  small  increments  in  u),  it  can  be  employed  advantageously  to 
provide  a first  approximation  to  the  source  location.  To  obtain  a more 
accurate  fix  on  the  source  location,  the  two  equations,  Incorporated  in 
(4-1),  can  be  treated  separately,  and  independently  solved  for  9 as  a 
function  of  0 based  upon  some  initial  estimate  of  0e . Using  an  initial 
estimate  of  0„  (say  0*),  then  two  independent  values  of  9 will  be  obtained 
trom  (4-1).  By  observing  the  two  values  for  9,  and  noting  which  value  is 
the  larger,  a second  estimate  for  0„  may  be  made.  Again  two  values  for  9 
will  be  obtained  and  the^e  can  be  compared  with  each  other,  and  with  the 
earlier  results  obtained  for  0: , to  arrive  at  a still  more  refined  estimate 
of  0O . In  this  iterative  manner,  the  solution  for  90  and  0t  can  be  con- 
verged upon  in  a relatively  short  time  with  the  aid  of  high  speed  digital 
computers . 

Although  one  might  solve  equation  (4-1)  explicitly  for  0S , this  process 
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would  result  in  an  extremely  cumbersome  relation  for  all  except  the  most 
simple  geometry.  Therefore,  this  procedure  is  not  recommended.  What  is 
recommended  is  to  consider  the  relation,  given  in  (4-1),  as  two  independent 
equations  for  8 as  a function  of  0,  and  solve  these  equations  for  various 
estimates  of  , until  the  resulting  8 solutions  converge  to  the  sate  value 
(viz.  8C ) . An  algorithm  for  rapidly  converging  to  the  desired  solution 
(8„,0O)  will  now  be  given. 


4 . 3 Iteration  Algorithm 

To  begin  the  iterative  solution,  we  shall  first  require  an  esti- 
mate of  the  value  0,,  . This  estimate  reed  not  be  too  accurate  and  may  be 
obtained  from  a rather  crude  graphical  plot  of  the  solution  as  described 
earlier.  One  should  then  choose  two  values  of  0 equally  spaced  about  the 
best  estimate  of  0O . For  convenience,  we  shall  call  these  two  choices, 

0j  and  0_.  Using  these  values  for  0,  we  then  compute  8j,j  and  9al  (for 
0=0^  and  0^  and  0a2  (for  0«=0?),  using  the  two  equations  giver,  in  (4-1). 
Based  on  the  results  for  the  first  two  choices  for  0 (viz.  0 j and  0P),  we 
•hen  compute  the  best  estimate  for  0„  as  0s  where, 


03  “ 


= 0 


? 


( s " ^1.  a )#i  - ( 8P,  j - 0j , j )08 

( 0?.  ? ‘ , ? ) ' ( 0s.  i " ®l,  l ) 

(0?,p-pi.  *H*,-*Q 

(0?.  r'01,  a )'  ( 0?,  l " 01,  1 ) 


(4-2) 


Using  the  estimate  03  in  (4-1),  then,  the  values  for  05i  3 and  9Sf  3 
are  computed,  and  a further  estimate  for  0„  is  obtained.  This  iterative 
process  is  continued;  whereby,  newer  and  more  refined  estimates  for  0C  are 
obtained  from  the  general  formula, 


(0?,  ,-6i,  .)  <0. -*„-i ) 

' 0?.  n " ®1,  n ) " ( 0&  b-  1 “01,  n-  1 ) 
n-1  n 02, 

= 0,-K-l)  (0?-0i)  n (4-3) 

(0?.  t '01,  1 )"(0r-  l-i) 
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Subsequent  to  each  estimate  of  0O , one  may  test  the  result  against 
a predetermined  tolerance  A0  and  AO  as  follows, 


| 0,.  + l - 0„  I < A 0 

I §2,  ti  + 1 “ 6j,  r + l ! < A 0 


(4-4) 


When  the  above  criterion  is  met,  the  process  may  be  terminated  and 
the  value  of  (0O,0O)  selected  as, 


(0O,0.)  - (0n+1  A + j)  (4-b) 

That  such  a solution  will  converge  in  the  limit  as  n—  ® should  be 
evident  from  the  general  geometry  and  from  (4-3).  For  most  practical 
applications  the  convergence  should  be  quite  rapid,  requiring  only  several 
iterations  (dependent,  of  course,  on  the  tolerence  values  A9  and  A0  and 
the  problem  geometry). 


Ill 

SOURCE  DYNAMICS  DETERMINATION 
OVER  SPHERICAL  SURFACES 


It  is  well  known  that  the  first  order  effect  of  relative  motion,  on 
the  temporal  structure  of  a signal  being  propagated  through  a transmission 
medium,  is  manifested  as  a time  scale-factor  change  on  the  signal  being 
transmitted.  That  is,  the  signal  is  either  compressed  rr  expanded  in  time 
(usually  by  a very  small  fraction)  as  a result  of  the  relative  motion 
between  the  source  and  the  receiving  sensor.  Receiving  sensors  in  the 
transmission  medium  will,  then,  measure  the  compressed  (or  expanded)  signal 
instLau  of  the  true  signal  which  was  transmitted.  And  the  signal  measured 
at  each  receiving  station  will  vary,  depending  on  the  velocity  of  the  source 
and  the  angle  between  the  velocity  vector  and  the  transmission  vector  to 
the  receiving  sensor,  liy  measuring  ar.d  noting  he  difference  in  time  scale- 
factor  shift  (dopplor  shift  in  the  case  of  C'W  signals)  between  two  (or  more) 
pairs  of  receiving  stations,  it  is  then  possible  to  uniquely  determine  the 
dynamic-'  (course  and  speed)  of  the  source. 

The  purpose  of  the  present  section  will  be  to  develop  the  algorithms 
which  express  the  source  dynamics  (course  and  speed)  in  terms  of  the  time 
scale-factot  shifts  which  are  measured  between  two  pairs  of  recelvii % 
sensors.  The  geometry  hypothesized  will  be  a spherical  surface  (presumably 
the  surface  of  the  earth)  on  which  the  source  and  receiving  stations  are 
located . 

1 • 0 PROBLEM  FORMULATION  FOR  A RECEIVING  SENSOR  PAIR 

Figure  r>  depicts  the  relative  geometry  of  the  s-'irce  dynamics  vector 
and  the  propagation  vectors  to  a pair  of  receiving  sensors  located  on  the 
surface  of  a sphere.  The  locations  of  the  source  and  receiving  sensor 
stations  are  presumed  known  and  the  signal  transmission  paths  are  assumed 
to  be  great  circle  routes. 

The  time  scale-factor  shit,  % , realized  at  the  receiving  sensor 
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TRUE  NORTH 


FIG.  5 VECTOR  GEOMETRY  OF  THE  PROBLEM 


(due  to  source  motion)  will  be,* 


where, 


(1-D 


v 1 s the  speed  of  the  so>  in  the  transmission  medium 
8 is  the  course  of  the  source  relative  to  true  North 
Bo,  is  the  direction  of  the  great  circle  path  between  the 

source  and  the  xc^  receiving  sensor  (»-elative  to  true  North) 
c is  the  velocity  of  propagation  in  the  i >dium 
(Th  ' above  approximate  relation  will  be  highly  accurate  providing  that 
v « c . ) 


* see,  for  example,  A.  A.  Gerlach  "Theory  and  Applications  of  Statistical 
Wave-Period  Processing",  Gordon  and  Breach  Science  Publishers  Inc., 

Vol . I,  Chap.  3,  p.  98;  1970.  ' 
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The  difference  in  time  scale-factor  shift  that  will  be  measured 
between  any  two  receiving  sensor  stations  y and  x,  will  therefore  be. 


X ( V , g ; 8 o T , g g ..  ) - f'y  - 6, 

^ - fcos(6„y-0)-cos(8o,-B)] 
c 

- 2 ^ sin^(g„,-goy)sinf^(B0T+eo/)-g]  (1-2) 


Tlie  above  relation  contains  two  unknowns,  v and  0,  and,  therefore, 
cannot  be  employed  to  uniquely  determine  the  source  dynamics.  However, 
the  one  unknown  can  be  explicitly  solved  for  in  terms  of  the  other  to 
determine  the  relationship  between  the  two.  Thus,  we  may  obtain  from  (1-2), 


and 


e 5yx 

2 sin  -60y  )sin[^(PriI+B0y  )-g] 


0 = ^(0ox+0oy )-sin-1{ 


c 6. 


y * 


\ 

2vsin^(B0X -B„y ) J 


(1-3) 


(1-4) 


As  a consequence,  if  one  or  the  other  of  the  unknowns  can  be  reasonably 
well  estimated,  the  other  can  be  computed  from  either  (1-3)  or  (1-4). 
However,  it  will  be  observed  that  bivalued  ambiguity  does  exist  (in  the  case 
of  p ,1  on  account  of  the  anti- trigonometric  relation.  In  certain  applica- 
tions this  ambiguity  might  be  removed  by  physical  considerations.  In  any 
event,  two  (or  more)  pairs  of  receiving  sensors  will  be  required  to  uniquely 
solve  for  the  source  dynamics. 


2.0  DETERMINATION  OF  THE  PROPAGATION  VECTORS' 

before  proceeding  w,th  the  unique  solution  of  the  source  dynamics,  it 
will  be  well  to  develop  tha  relations  for  the  great  circle  propagation 
vectors,  g„„,  between  two  points  ( 60  , > and  (d„,',\). 


To  determine  the  angle  of  propagation  between  the  two  relevant  points 
on  the  surface  cf  the  sphere  we  will  first  compute  the  slope  of  the  great 
circle  path  between  the  two  points.  From  equation  (2-31 ) in  Section  II, 


we  may  write  the  equation  for  the  great  circle  curve  through  the  two 
points  (6.  A)  and  ( 0„  ,0X ) as, 

q _ tlin-i  jtan6°sin(fl-0, )+ tan 8X sin (0,-0)  j 
U"tan  | sin  (0e  A)  j 

where  the  angles  80 , 8,  and  0O  ,0,  are  given  in  North  Latitude  and  West 
Longitude  coordinates  respectively. 


Now,  the  slope  of  this  curve  may  be  readily  shown  to  be, 


tan[a(0;  0„  ,0„  ,6,  ,0,)] 


i de 

cos6  d0 


tan8cco8(0-0, ) - tan8, 008(0.  -0) 

^sin v(ifie  -0, )+  [tan60  ain(0-0,  )+tan6,  sin(0e  -0)]  ~ 

Therefore,  the  angle  a,  of  the  propagation  vector  at  the  source  location 
(relative  to  the  Latitude  parallel),  will  be, 


»<*.) 


tan 


(sln8n  cos  (0e  -0, ) -cos8n  tan 8, 
| sin(0o-0, ) 


(2-3) 


Since,  now,  a will  range  over  ± 90  degrees  we  shall  need  to  Incorporate 
the  following  additional  algorithm  for  the  computation  of  the  propagation 
angle  B„,  (relative  to  true  North), 


i o,(0o  )+90°  when  0,  0„  i 

Be,  J (2-4  > 

(r,(0o)+27 0°  when  0,^  0,  j 


It  will  be  of  interest  to  compute  the  direction  of  the  propagation 
vector,  8„„,  from  the  xttl  receiving  station  toward  the  source.  Kr.  i (2-2) 
we  have, 


alc>  ) „ tan-i  J CP.?A  t an 0,_-_s  In  5, cos(g,  -0,  } { 
(v,)  C3n  \ sin(0o -0, ) j 

and,  the  angle  g„„  (relative  to  true  North)  becomes, 

| '*(?,)  + 270°  when  0,  0o 


(2-5) 


Bxc 


"*(01  ) +■  90°  when  0O  e 


J 


(2-6) 
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For  the  particular  situation  when  0,  = 0O , the  following  algorithms 
apply, 

i 180°  when  0,  < 8„  j 


®o  * = 


and, 


e 


0°  when  8e  <*  8, 

0°  whe  0,  < 0C 

180°  when  8e  < 8, 


(2-7) 


(2-8) 


It  is  significant  to  note  that,  in  general,  the  angle  B„0  is  not  simply 
equal  to  the  angle  B0«  minus  180°.  As  a matter-of-fact,  under  special 
circumstances,  B«0  will  be  equal  to  Bo* • This  apparent  dilemma  is  due  to 
the  nature  ot  the  great  circle  arc  over  the  surface  of  the  sphere. 

Some  symmetries,  relative  to  the  above  propagation  vectors,  are  as 
follows:' 


From  (2-3),  it  may  readily  be  shown  that, 

o {(180+$,  ) - (0O  -0,  1 ; - 0,  , (180+0,  ) - (0O  -0,  ) , 0,  ,0, )} 


“ 3(00 ; 8o  >0o  1 0«  ,0* ) (2-9) 

fH0,  - (0o  -0, ) ; 60 ,0,  - (0o  -0, ) , 9,  ,0, )} 

- -o(0o ; 0o , 0.,  9,  A)  (2-10) 

tr,  0B  ; 0o,0o -6,,0,  + 180 J 

- 3(0o;0o,0c ,0. ,0.)  (2-11) 

and , 

3‘0o  ; 00  ,00 , 0,  ,0«  +2  (0O  -0,  )J 

= -<*(0o  ;6o  ,00 ,0*  ,0, ) (2-12) 
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And  from  (2-5),  it  can  be  shown  that, 


; - e0,(  18O+0, )-  (0O  -0, ) , e,  ,0, } 


= -*(0*;eo,0o,e,,0x)  (2-13) 

^0* ; 00 ,0*  - (00  -0* ) . 6*  ,0x  j 

- -«(0,;eo,08,e,,0,)  (2-i4) 

^0,+18O;0o  ,0O  ,-8,  ,0,  + 18o!- 

j, 

= -G(0,,8o,0e,O,,0,)  (2-15) 

fyJ0,+2(0o  -0,  ) ; 0O  ,0O  , 6,  ,0,+2(0o  -0,  )| 


" -^(0*;6o»0o.0*»0*)  (2-16) 

Equations  (2-9)  and  (2-11)  demonstrate  the  convergence  zone  property 
of  the  spherical  geometry  for  a source  located  in  the  extreme  far  field. 

Thus,  the  apparent  sensor  pair  aperture  angle  for  a source  located  in  the 
extreme  far  field  (on  the  opposite  side  of  the  hemisphere  from  the  sensor 
locations)  will  be  the  same  as  for  a source  located  in  the  corresponding 
image  position  on  the  near  side  of  the  hemisphere.  This  convergent  zone 
property  will  permit  the  source  dynamics  measures,  of  a source  located  in 
the  extreme  far  field,  to  be  resolved  with  the  same  precision  as  would  be 
realized  in  the  corresponding  image  location  in  the  near  hemisphere.  More 
will  be  said  about  this  feature  of  the  spherical  geometry  in  Section  V. 

3.0  SOURCE  DVNAMICS  ALGORITHMS 

In  the  earlier  analyses  we  determined  that  the  computation  for  the 
source  dynamics  (course  and  speed)  would  require  the  solution  of  two  sets 
of  nonlinear  simultaneous  equations  of  the  form  given  in  (1-2).  Thus,  if 
we  assume  that  two  pairs  of  receiving  sensors  were  receiving  and  simul- 
taneously measuring  the  time  scale-factor  shift  from  the  transmitted  source 
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signal,  the  two  resulting  nonlinear  equations  involving  the  unknown 
dynamic  variables,  v and  0,  would  be, 


^y*  - 2 ~ sln^j ( 80  , -8 0 y ) sin  [ ^ (g 0 „ +0 6 y ) - 8 ] 

\,t  = 2 - s in% (B o i ■8o*)sin[^(beI+g0»)-g] 
c 


(3-1) 


where  the  two  pairs  of  receiving  sensors  are  at  the  known  locations  (0«,0, ), 
(0y,?y)  rand  (0t  ,0t ) * (0*  »0»  ) i and  the  corresponding  propagation  vector  angles 
are  8mSoy  and  floi»8o«»  respectively. 


For  convenience,  we  shall 


let, 

sin^(Boy  -8o») 
sinlj(eow  -0O, ) 


and, 

^ * ).  (®o  x +8  o y 't'Bo  z +8  o v ) 

n “ Ueei+6o»'BeK-Poy) 
From  (3-1),  (3-2),  and  (3-3)  then, 


(3-2) 


(3-3) 


v sinfci-B-'H)  » -c6y,/2byl 
v sin(:-B+n)  = -c6,,/2b>. 


(3-4) 


Making  use  of  the  identity,  sin(a±b)  = sina  cosb  ± cosa  sinb,  the  above 
expression  reduces  to,  simply, 


v sin(J-g) 
v cos('-8  ) 


by  x % t+by  , f>rx 
by  i b»  r cost] 


by  ,b,  i sinrj 


(3-5) 
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The  above  two  nonlinear  equations  in  the  unknown  variables,  v and  g, 
may  be  readily  solved  to  give, 


8 


•tan 


.1 1 % t+b|rt,5yt 

*•  by  X 5«  l -h*  J by  X 


tanr) 


(3-6) 


and 


c 1 (bY  x &HJ.)  + (!?¥  i.3t)  ‘'"foxakx  i by  i Pi  i C08^rll, 

2by,bwt  sin2T] 


(3-7) 


Thus,  the  source  dynamics  (course  and  speed)  may  be  explicitly  deter- 
mined from  the  measured  time  scale-factor  shifts,  fy*  and  5*,,  and  a knowledge 
of  the  locations  of  the  source  and  the  receiving  sensor  stations.  To 
avoid  ambiguity,  which  may  result  from  the  "sign"  of  the  trigonometric 
functions,  some  care  should  be  exercised  in  the  use  of  (3-6)  and  (3-7). 


It  will  be  noted  in  (3-7)  that  v will  always  be  a real  number;  however, 
it  is  possible  for  v to  take  on  negative  values.  In  the  event  that  v is 
negative,  the  simple  interpretation  is  to  take  onl”  the  absolute  value  for 
v and  add  ± 180  degrees  to  the  value  of  8 to  make  it  fall  within  the  range 
0 < 0 360.  To  avoid  the  ambiguity  of  a negative  speed,  a simple  pre- 

caution is  to  choose  the  sensor  nomenclature;  such  that,  by,,b,,  and  sin2r] 
are  all  positive  or  that  two  of  these  terms  are  negative.  For  many  appli- 
cations this  can  be  accomplished  by  insuring  that 


providing  that, 


8C  x ‘ Boy  Bo  I Sol 


Sol  "F  Sow  ~ Bo*  * Boy  360 


(3-8) 

(3-9) 


Since  three  receiving  sensor  stations  are  sufficient  to  produce  two  pairs 
of  receiving  sensors,  we  can  let  0,  = 9y  and  0t  = 0y  . In  this  event, 


'or 


and  (3-3)  becomes, 

i 

I 8 = \ (Bo  x F 280y  +60W  ) 

I 71  = i(8.»-6ox)  • 

In  all  other  respects  the  relations  are  unchanged. 


(3-10) 


(3-11) 
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IV 


RESOLUTION  IN 
SOURCE  LOCALIZATION 

In  Section  II  we  developed  the  algorithms  (or  analytical  relations) 
for  computing  the  location  of  a signal  source,  based  on  the  time  register 
measure  between  two  pairs  of  receiving  sensors  whose  locations  were  known. 
Hecore  the  practicality  of  these  algorithms  can  be  properly  assessed, 
however,  it  will  be  necessary  to  determine  the  resolution  of  this  locali- 
zation in  terms  of  the  resolution  of  the  basic  time  register  measure  T,  In 
other  words,  we  wish  to  ascertain  the  sensitivity  of  the  localization 
incisure  to  the  inherent  resolution  of  the  measure  in  time-difference- in- 
arrival of  the  source  signal  at  the  receiving  sensor  pairs.  When  this  is 
accomplished,  we  can  suitably  outline  an  area  on  the  surface  of  the  sphere 
in  which  we  can  be  assured  that  the  source  is  known  to  exist  within  a 
definable  degree  of  confidence.  This  we  shall  now  proceed  to  do. 

1 • 0 PROBLEM  FORMULATION  FOR  A RECEIVING  SENSOR  PAIR 

The  approach  to  the  problem  will  be  one  of  determining  the  sensitivity 
cf  the  equi- time-register  curves  to  the  nondimensional  parameter  u.  Since 
the  parameter  u is  directly  related  to  the  time-register  measure  T,  as  given 
in  H-4)  of  Section  II,  we  shall  then  be  in  a position  to  determine  the 
sensitivity  of  localization  to  the  measure  T.  Our  plan  will  be  to  ascertain 
both  the  translation  distance,  ds  influenced  by  an  error,  du , in  the  measure 
u,  and  tlie  direction  of  this  translation.  When  this  is  accomplished  for  two 
pairs  of  receiving  sensors,  we  shall  then  be  in  a position  to  outline  an 
area  on  the  surface  of  the  sphere;  in  which,  the  source  will  be  localized 
for  a specified  error  in  time-register  measure. 

Since  the  geometry  of  the  problem  is  independent  of  the  coordinate 
frame  of  reference  of  the  sphere,  we  shall  conveniently  choose  the  reference 
frame  given  in  Section  II,  3.1. 

1 . 1 Fundamental  Equations 

The  fundamental  equations  which  we  wish  to  solve  may  readily  be 
shown  to  be, 
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ds  _ p ! d6/du  1 

du  Jl  + (d9/d0)a/cos'^ 


end , 


8u 


d0/du 

tan 

t d0/du  j 


fc£s9_\ 

-d6/d0J 


(1-2) 


where  ds  Is  a measure  of  distance  over  the  surface  of  the  sphere  and  Bu  is 
the  angle  along  which  the  distance,  ds,  is  to  be  taken.  Since  the  angle 
8U  can  be  rather  ambiguous  due  to  the  inherent  nature  of  the  transcendental 
functions,  it  will  be  best  to  determine  this  angle  from  a knowledge  of  the 
geometry  of  any  specific  problem.  In  essence,  it  is  the  direction  normal 
to  the  tangent  of  the  equl-tlme-register  curve  (at  the  point  in  question), 
and  toward  the  curves  of  higher  u.  Thus,  the  measures  ds  and  gu  define  an 
error  vector  on  the  surface  of  the  sphere;  where  ds  is  the  scalar  error 
magnitude  and  gu  is  the  direction  of  the  error  vector.  It  should  be  evident 
that  the  error  magnitude  will  be  simply, 


ds 


du  ds  . 

dl  d^  dT 


(1-3) 


It  will  be  noted  that, 


And  for  the  geometry  given  in  Section  II,  3,1,  Bu  will  be, 


(1-4) 


(1-5) 


where  Bu  is  expressed  in  angular  degrees  relative  to  true  North.  Now, 
the  slope  of  bu»  as  given  in  (1-5),  is  given  in  (3-14)  of  Section  II. 
Consequently,  it  will  be  convenient  to  utilize  the  abo''“  expression  for 
Bu  in  the  analysis  to  follow;  in  lieu  of  the  expression  given  in  (1-2). 
Since  it  will  be  well  to  examine  the  geometry  of  any  specific  problem;  :o 
determine  the  true  direction  of  bu«  the  actual  form  of  this  expression  in 
the  subject  analysis  is  somewhat  arbitrary;  as  long  us,  the  final  results 
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are  not  misapplied.  Consequently,  for  convenience,  wf  shall  use  the  form 
given  in  (1-5)  in  the  analysis  to  fellow.  In  general,  this  will  not  lead 
to  ambiguous  results;  as  long  as,  care  is  taken  in  the  employment  of  the 
results  to  any  physical  problem. 

1. 2 Graphical  Interpretation  of  the  Error  Vector 

At  this  point  one  may  obtain  a general  feel  for  the  quantity  of 
the  error  vector  through  a study  of  the  curves  shown  in  Figures  2,  3,  and 
4 of  Section  II.  For  any  given  point  on  the  family  of  curves,  one  may 
draw  a vector  from  this  point  normal  to  the  curve  in  the  direction  of 
increasing  u.  The  distance  between  two  adjacent  curves  divided  by  the 
incremental  difference  in  the  values  of  u will  be  the  approximate  magni- 
tude of  ds/du.  And  for  the  subject  geometry,  the  value  of  du/dT  is  simply, 


dii 

dT 


d-6) 


where  D is  the  great  circle  distance  between  the  two  receiving  sensor 
stations  and  c is  the  velocity  of  signal  propagation  in  the  transmission 
medium. 


2.0  SOLUTION  OF  THE  BASIC  ERROR  EQUATIONS 
2 . 1 General  Solution 

Making  use  of  equations  (3-4)  and  (3-14)  In  Section  II,  3.1,  we 
can  solve  the  fundamental  equations  given  in  Section  1.1  above  to  yield 
the  following  relations. 

. n [sinr0,  cos”"  ,i  0,  - (sin20x  - sin4u0x  )cos20x  cosr0]  0„ 

dj  „ __ --  - — - — - ■ - (2-1) 

u (3in20,  cos'  u0* -cos30,  sin‘ u0,  cos20)  (sin30x  -sin2u0x ) (cos3u6x  -cosa0x  cos20)k 


(d0/d0)  5 _ f sin"  9.  cos4  '~i~,  - ( sin3  6,~ 
cos' 0 V (cos'3** O, -cos'  0j cos“0) 


■sin4  u9,  )cos50,  cos20 


(sin30xcos''u0x  -cos20«  sin"u0jt  cos30) 


(2-2) 
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ds 

du 


a 


P0x 


slng6xcos‘tu8x  -(sin58x  -sln*u8x  )coss8x cos*^ 

(slng0, -sin2u  8X ) (sins8, cosSu8j -coss0x  sin3a  8X  cos ?0 


(2-3) 


= tan 


.1 


cos?0x  sinu 0, sin0cos0 


(2-4) 


Although  the  above  expressions  3re  rather  complex,  they  do  simplify 
considerably  for  particular  values  of  the  variables  u and  0.  For  example, 
when  u is  equal  to  zero,  the  two  relevant  relations  become, 


ds 

duJ 


/ 


1 4-cot3 0X  sin 20  p6x 


u" 


(2-5) 


Puj  - o 

u=0 


For  0 equal  to  zero  and  ninety  degrees,  the  expressions  become. 


ds 

du 


0=0 


pex 


(2-6) 


■»] 


6u 

0=0 


and 


P0X  COSii0x 

ds  1 _ — — — 

du  0=n/2  " V sln<,ex-sinyU0x 


(2-7) 


0=tt/2 


In  all  of  the  above  situations,  the  vector  direction,  Bu>  is  equal  t0 
zero  relative  to  true.  North.  That  is,  the  slope  of  the  equi-time- register 
curves  is  zero  at  the  relevant  points. 
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2.2  Solution  for  Small  0, 

Of  particular  interest  will  be  the  case  where  the  angle  6,  is 
sufficiently  small,  so  that,  the  sine  of  the  angle  can  be  replaced  by  the 
angle  (expressed  in  radians)  without  undue  error.  This  is  not  too  great  a 
restriction,  since  this  criterion  can  be  met  for  sensor  pair  spacings  which 
are  equal  to  or  less  than  about  1,200  nautical  miles  over  the  earth's 
surface.  With  this  assumption,  then,  equations  (2-3)  and  (2-4)  reduce  to, 


and , 


ds 

du 


(i-us)s9,2+(i-e!ts)(i-u4  e,a)sinV 

( 1-u"  ) [ l-ua  (coss0+0, 3sin20)  ] 


(2-8) 


Bu  = tan 


-1 


u(l-G, H)sin0cos0 


\f  ( l-u.3(cos204-9x ssin20) ] [sin5<H0*S(cosS0-u2)] 


(2-9) 


For  convenience,  we  can  now  study  the  problem  in  the  near  field 
(including  the  extremely  far  convergent  zone)  and  the  relatively  far  field 
(remote  from  the  receiving  sensors). 


2.2.1  Near  (and  Convergent  Zone)  Field  Solution:.  In  the  near 
field  and  in  the  extremely  far  field  (or  convergent  zone-',  where  the  sine 
of  the  angle  <t>  can  be  replaced  by  the  angle  itself  (or  TT  minus  the  angle), 
the  relevant  equations  approximate, 

~ * p y e,3+0a/(l-ua)a  (2-10) 


Bu  « tan-1 


u0 

yaVW+Va-uV 


(2-11) 


(NOTH:  For  the  convergent  zone  area,  the  variable  0 should  be  replaced  by 

T-0.) 

2.2.2  Far  Field  Solution;  For  the  far  field  solution  ^remote  from 
the  sensor  locations),  the  relevant  relations  may  be  approximated  as, 
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(2-12) 


, p sin  0 

ds  p 1 ==  , _ ; 

du  “ y (l-u2)(l-uscos20) 

« tan-1 

3.0  TRANSFORMATION  OF  THE  COORDINATE  FRAME  OF  REFERENCE 

In  the  earlier  analysis  of  the  basic  error  relations,  we  developed  the 
equations  for  a particular  geometry  of  the  receiving  sensors  in  the  normal 
spherical  coordinates.  Since  we  wish  to  apply  these  results  to  any  geometry 
over  the  surface  of  the  sphere,  we  shall  now  transform  the  specific  spherical 
coordinate  measures  to  a new  measure  space  which  is  directly  related  to  the 
physical  geometry  of  the  problem  and  independent  of  the  frame  of  reference 
of  the  sphere. 

To  accomplish  the  above  objective,  the  new  coordinate  frame  of  reference 
will  b ■>  vhe  midpoint  of  the  great  circle  curve  which  joins  the  two  receiving 
senso-  i.  (In  the  above  example,  this  point  is  the  point  ft*O,0«O.)  We  then 
propose  to  employ  two  new  variables,  R and  a,  where  R is  the  great  circle 
range  from  the  reference  point  to  the  source  position,  and  a is  the  angle 
between  the  great  circle  path  R and  the  normal  to  the  great  circle  path 
between  the  two  receiving  sensors;'  taken  at  the  reference  point. 

3 . 1 Transformation  Equations 

For  the  example  problem,  the  normal  to  the  great  circle  path 
between  the  two  receiving  sensors  at  the  reference  point  is  simply  the 
equator,  whose  slope  is  zero.  The  slope  of  the  great  circle  path  R (at  the 
reference  point)  may  be  obtained  directly  from  (2-3),  Section  III.  Thus, 
we  have  the  relation 


u cos  0 

i/  1*U2C0S  20 


(2-13) 


tan  a *=  tan0/sin0  (3-1) 

The  great  circle  range  R may  be  obtained  directly  from  (1-5)  of 
Section  II,  to  give  the  relation, 

sin(K/2p)  = ^ (1-COS0COS0)'  (3-2) 
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From  the  above  two  nonlinear  equations  in  the  new  variables  R and  a, 
we  may  solve  explicitly  for  the  new  variables.  The  useful  interrelations 
between  the  four  variables  have  been  derived  and,  for  convenience,  they 
are  displayed  in  Table  1 . 


TABLE  1 INTERRELATIONSHIP  BETWEEN  VARIABLES 


sin'  C 

sr 
cos  c 

tan3C 

R 

P 

l-cos2Osin30 

cos30sin30 

tans8  + slny0 
cos20 

a 

tan2  8 

sin20 

tan30 

tan20Fsin30 

tan30+sin30 

sin30 

0 

sin?(R/ p) sin3# 

l-sin3(R/p)sin3Q' 

sin3(R/p)sin3Q' 

l-sin3(R/p)sinJQ' 

0 

sinp(R/p)cos?Q' 

cos2 (R/ p) 

tan3  (R/p)co33» 

l-sin3(R/p)sina# 

l-sin3(R/p)sin2# 

3. 2 Solution  for  the  Parameter  u 

With  the  introduction  of  the  two  new  variables,  the  parameter  u 
is  no  longer  independent,  but  implicitly  defined  by  these  variables.  The 
solution  for  u in  terms  of  the  new  variables  is,  however,  not  a simple 
function,  but  fortunately  can  be  derived  explicitly. 

From  (3-2),  Section  II,  we  may  write, 

X4-(sin29I+sin*’t3fcoss8,cosP0sin?0)XE+sins9Isinc'8  = 0 (3-3) 

And  using  the  relationships  given  in  Table  2 we  obtain, 

X4-[sln3  ^ + sin20„ (1-sin2  ^ cos2^) ] X2+sin3 8,  sin2  ^ sin2a  = 0 (3-4) 

where,  X2  » sin2u8,  5 sin3*),  (3-5) 

Equation  (3-4)  has  the  solution, 
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X2  - % [b-  /IP -4a  ] 

where, 

b * sin*  — + sin*^  (1-aii*2  — cos2^) 
P P 

a =»  9in20x8in2  — sin3** 


[NOTE:  Only  the  smaller  of  the  two  solutions  for  X3  need  be  considered  in 

(3-6),  due  to  the  restriction  imposed  by  (3-5).] 

When  0*  is  sufficiently  small  so  that  sin6„  can  be  approximated  as 
0* , equation  (3-4)  becomes, 

e,V-[sin2  | + 0X 2 (1-sin 3 * cosaGr)Jua+sina  ~ sina<*  - 0 (3-8) 

This  latter  equation  will  be  solved  for  boch  the  near  (including 
convergent  zone)  field  and  the  far  field. 

3.2.1  Near  (and  Convergent  Zone)  Field  Solution:  For  the  situation 

where  R/p  (or  rr-R/p)  is  sufficiently  small  so  that  oina(R/p)cosaa'  « 1 , 
equation  (3-8)  may  be  approximated  as, 


where, 


u4-(Sa+i)u3  + §a«inac*  - 0 


5 ■ sin(R/p)/9i  ■ 2p  sin(R/p)/D 


(3-10) 


The  solution  to  (3-9)  is, 


/ - - J l+?4+2§3cos2<y  ] 


(3-11) 


3.2.2  Far  Field  Solution:  When  sin(R/p)  is  sufficiently  large 


so  that, 


sin* (R/p)»  Da/',4pa+Dacobacv) 


(3-12) 


equation  (3-8)  may  be  approximated  as, 
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u‘;-rus+e®sinaa  * 0 


(3-  13) 


The  solution  in  this  case  is 

us  - \ f[l-J  l-(2sinc*/!()3' j 

• a r t i 9 in  a . sin  o'  , i , » . i ^ 

= sLn  tv  [1+  — — + ? — — + ...]  (3-14) 

^ -* 

4 . 0 SOLUTION  OF  THE  ERROR  VECTOR  FOR  A RECBxv’ING  SENSOR  PAIR 

Using  the  relations  given  in  Table  2,  and  '.he  solutions  for  u given 
in  the  last  section,  we  may  readilj  compuce  the  error  vector  (for  a given 
receiving  sensor  pair)  in  terms  of  the  new  variables  R and  tv.  This  will 
be  done  for  both  the  near  and  far  fields. 


4.1  Near  (and  Convergent  Zone)  Field  Solution 

For  the  near  field  solution  we  shall  assume  that  sin3  - « 1 and 

P 

employ  the  relations  given  in  (2-10),  (2-11),  and  (3-11)  along  with  the 
data  provided  in  Table  2,  to  obtain. 


and, 


ds 

du 


V 


1+^+2^00320  +(1-^)7  1+5* +2?,3co8  2c* 
l+ff4-2f'"ain3o  +(1-=(V  l+e;4+2f,3cos2o 


(4-1) 


tan  <?u  *■ 


± 


?3cos3o  [ 1+^3-  J 1+^+2^008  20  I 

l+^sin3^  +!r3(l-3sin3o)+(l-ff3sin2a)^  I+^+IFcoTTtT 


(4-2) 


(NOTE:;  Due  to  the  change  of  variables,  the  angle  specified  as  no  longer 
will  give  the  direction  of  the  error  vector  relative  to  a true  North  refer- 
ence. A new  variable  is  therefore  employed  in  place  of  the  former 
variable  gu  which  was  applicable  only  to  the  earlier  specific  example.  In 
the  present  analysis,  the  angle  »u  will  be  orthogonal  to  the  measure  of  the 
error  vector  direction.  The  ingle  is  measured  relative  to  the  slope  of  the 
great  circle  curve  joining  the  point  (R,aO  and  the  midpoint  between  the  two 
receiving  sensors;  the  referenced  slope  to  be  taken  at  the  point  (R,a). 

The  geometry  will  be  demonstrated  later  in  Figure  10.) 
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The  above  equations  can  be  greatly  simplified  by  considering  limited 
ranges  for  the  value  of  ?. 

4.1.1  Case  Where  2 £ * : For  the  situation  where  the  variable  ? 

is  equal  to  or  greater  than  two,  equations  (4-1)  and  (4-2)  reduce  to, 


Making  use  of  (1-6),  equation  (4-3)  becomes, 

~ yl+[Dcoflo/2p8in(R/p))8'  (4-5) 

dT  Dcoo<v 

It  should  also  be  noted  that,  when  * is  large,  the  direction  of  the 
error  vector  is  approximately  orthogonal  to  the  angle  o;  since,  under  the 
circumstances,  ! ou  ! « j n 1 . 

4.1.2  Case  Where  ^ - 1:;  For  the  situation  where  the  variable 
f Is  equal  to  unity,  equations  (4-1)  and  (4-2)  become  simply. 


and, 


^ - D/  /T 
du 

«U  “ 0,/2 


(4-6) 

(4-7) 


Thus,  at  the  particular  range  R=l)/2,  the  scalar  magnitude  of  the 
error  vector  is  independent  of  the  angle  a.  In  addition,  the  direction  of 
the  erro’’  vector  is  orthogonal  to  the  angle  or/2. 


4.1.3  Case  Where  g s 1/2:  When  the  variable  p2  is  small  compared 

with  unity,  equations  (4-1)  and  (4-2)  reduce  to, 
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u 2 y 


^ + ___i_££s_£ 

l^^sin2*  + ?*  (l-cos4 a) 


and, 


au  « ± tan 


.1 


F sino  coso 


/ l-^acQ8ao' 
Vl+^a-Ssin2®)*! 


)+^sinso 


(4-8) 


(4-9) 


Thus,  In  the  limit  as  F-O,  the  magnitude  of  the  error  vector  reduces 
to  (c/2)dT  and  the  direction  follows  the  great  circle  path  between  the  two 
receiving  sensors. 


4.1.4  Graphical  Presentation  of  the  Error  Vector  Data:  Graphical 

plots  of  the  error  magnitude  and  the  angle  are  illustrated  in  Figures  6 
and  7 for  the  near  (and  convergent  zone)  field.  The  data  are  displayed  in 
convenient  format  for  data  abstraction,  once  the  values  of  range  and  rela- 
tive bearing  to  the  source  are  known. 


4.2  Far  Field  Solution 

For  the  far  field  solution  we  shall  assume  that  D2/4ps«sin2(R/p) . 
Employing  the  relations  given  in  (2-12),  (2-13),  and  (3-14)  along  with  the 
data  given  in  fable  2,  we  obtain 


1 


2pssina(R/p) 


y4p2sin?  (R/c)cossQ'-D2  ( 1+cos2 (R/p)  ] sin4  a 


and, 


au  « ± tan’1 


co s"  (R/p)  [4p  sin^R/pJ+D^sin'o)  sin3^ 


4 cS  h i n2 ( R/ p ) co  sS»- D2 co  s 2 ( R/ p ) a i n*  o 


(4-10) 


(4-11) 


The  above  equations  will  be  suitable  approximations  of  the  subject 
measures  for, 


I *>  ! tan’1 


£ 


sin£R/_E2  I 

° f 


(4-12) 


Ove:  the  indicated  range  of  o,  then,  the  error  vector  relations  may 
be  approximated  as 
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~ ££  sin(R/o) 
dT  D co  so* 


(4-13) 


o-u  k i tan*1  !cos(R/p)tanry  j (4-14) 

Graphical  plots  of  the  above  functions  are  illustrated  in  Figures  8 
and  9 for  convenience  in  data  abstraction. 

5.0  SOURCE  LOCALIZATION  ERROR  ALGORITHMS 

In  this  section  we  shall  describe  the  general  procedures;  whereby,  one 
may  make  use  of  the  formulas  and  graphical  illustrations  developed  in  the 
earlier  sections,  to  rather  accurately  define  the  area  In  which  a source 
has  been  localized.  This  area-of-localization  concept  (rather  than  a 
single  point)  is  based  on  the  fact  that  the  time  register  measure  T lias  a 
limit  of  resolution  which  will  be  ±At/2.  Since  the  measure  of  the  variable 
T is  not  precise,  there  will  be  some  ambiguity  about  the  exact  location  of 
the  source  which  is  directly  lelated  to  the  uncertainty  (or  resolution)  in 
the  measure  of  T.  Jt  is,  therefore,  significant  to  be  able  to  reasonably 
well  define  the  influence  of  the  t measurement  inaccuracy  on  the  locali- 
zation of  the  source. 


5 . 1 General  Description  of  the  System  Geometry 

To  define  the  source  localization  error,  it  will  be  necessary  to 
know  the  location  of  the  two  (or  more)  receiving  sensor  pairs  and  the 
general  location  of  the  source  relative  to  the  receiving  sensor  pairs.  In 
the  subject  problem,  the  precise  location  of  all  receiving  sensors  has  been 
presumed  to  be  known.  The  general  location  of  the  source  will  be  simply 
the  solution  of  the  equations  provided  in  Section  II.  This  general  location 
is  based  upon  the  actual  measures  of  t achieved  at  the  two  (or  more)  pairs 
of  receiving  stations.  (Or,  alternatively,  a general  location  can  be  simply 
hypothesized  for  purposes  of  error  analysis.) 

based  on  the  location  specifics  described  above,  one  may  then  construct 
a map  (or  diagram)  oi  the  general  geometry  as  illustrated  in  Figure  10.  The 
basic  geometric  orientation  is  illustrated  in  Figure  10(a),  where  the  point  0 
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depicts  the  position  of  the  source,  and  the  points  Aj  and  depict  the 
midpoint  of  the  two  sensor  pairs.  The  lines  OAj  (or  Rj)  and  OAs  (or  R^) 
depict  the  great  circle  paths  (and  ranges)  from  the  midpoints  of  the  sensor 
pairs  to  the  source.  The  lines  AjBj  and  AgBj  represent  the  great  circle 
arcs  through,  and  orthogonal  to,  the  midpoints  of  the  sensor  pairs.  The  two 
angles  and  az  represent  the  angles  to  the  source  from  the  respective 
receiving  sensor  pairs.  The  two  lines  Aj  ^ ' and  Ag  ^Bg  ' are  constructed  to 
pass  through  the  point  0 and  parallel  to  the  planes  of  the  arcs  A^  and 
AgB?  respectively.  The  angle  i is  the  angle  between  the  subject  two 
reference  axes. 

Attention  is  now  directed  to  Figure  10(b)  which  is  a blow-up  of  the 
region  about  the  point  0 with  the  t'«>  lines  Aj  'Bj  1 and  Ag  ' as  the 
appropriate  reference  axes.  However,  before  proceeding  further  with  the 
localization  geometry,  we  shall  first  neea  to  abstract  the  pertinent  error 
vector  data  from  the  appropriate  equations  or  graphs. 

5 . 2 Error  Vector  Data  Abstraction 

From  the  system  geometry,  depicted  in  Figure  10(a),  range  and 
angle  to  the  source  (R*,^  and  R^.o-j)  from  both  pairs  of  receiving  sensors 
may  be  determined  (or  computed).  Using  these  data,  along  with  the  specified 
time  register  measurement  resolution  (±At/2)(or,  alternatively,  ±0T) , one 
may  then  employ  the  earlier  equations  (or  graphs)  to  determine  the  magnitude 
of  the  respective  error  distances  (Asj  and  Asa),  and  the  relevant  angles 
(oul  and  rtU2)  for  each  sensor  pair.  These  data,  along  with  the  reference 
axes,  A^  ^ ' and  Ag  ' (illustrated  in  Figure  10)  may  then  be  employed  to 
construct  the  localization  error  parallelogram. 

(NOTE:  A suitable  measure  for  the  time  register  resolution  will  be  the 
standard  deviation,  vr,  of  the  time  register  measure  T,  Fither  ±CT  or  a 
specific  error  increment,  ±At/2,  can  be  employed  in  the  error  analysis.) 

5 . 3 Construction  of  the  Localization  Error  Parallelogram 
Referring  to  Figure  10(b),  the  abstracted  angles  cyui  and  fvua  are 

employed  to  construct  a set  of  error  axes  C1D1  and  CSD?.  (NOTE;  The 
direction  of  the  error  angles  may  readily  be  deduced,  from  the  general 
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FIG'.  10  DETERMINATION  OF  LOCALIZATION  ERROR  GEOMETRY 
FROM  SYSTEM  GEOMETRY  AND  ERROR  VECTORS 
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system  geometry,  by  noting  which  side  of  the  great  circle  path,  AB,  that  the 
source  is  located.)  Using  the  error  reference  axes,  then,  we  know  that  the 
direction  of  the  error  vector  is  orthogonal  to  these  axes.  As  a consequence, 
one  may  construct  the  error  parallelogram  by  drawing  appropriate  lines  which 
are  both;  parallel  to  the  error  axes,  and  spaced  a distance  ±As/2  from  the 
error  reference  axes.  The  two  sets  of  parallel  lines  are  extended  until 
they  intersect;  thus,  forming  the  error  parallelogram  illustrated  in  Fig.  10(b). 
(NOTH:'  The  commonly  referred  to  error  ellipse  may  be  constructed  within  the 
error  parallelogram  as  illustrated  in  the  figure.) 

Of  great  significance  is  the  angle  cp  of  intersection  between  the  two 
error  axes.  This  angle,  along  with  the  error  magnitudes  Asj  and  As2 , will 
determine  the  dimensions  and  geometry  of  the  localization  error  parallelo- 
gram. 


From  Figure  10(b)  one  may  see  that  the  angle  cp  is  simply, 

Cp  «=  f -Ouj-Ohs!  (4-15) 

The  lengths  of  the  two  sides  of  the  parallelogram  are  respectively, 

Ast/sincp  and  As2/slncp.  The  area  of  the  parallelogram,  thus,  becomes, 

Aren  - AsjAsg/sincp  (4-16) 

and  the  length  of  the  longest  diagonal  of  the  parallelogram  is, 

l _\/(Asi  )‘>(As, )a+2As1  Ass  coscp  (4-17) 

sincp 

The  importance  of  the  angle  cp  should  be  quite  evident.  It  can  be  sewn 
that,  the  resolution  of  localization  degrades  as  the  angle  cp  approaches 
zero  (or  tt)  , and  improves  as  cp  approaches  tr/2.  Optimum  localization 
resolution  occurs  when  the  source  is  located  along  the  two  great  circle 
arcs  AjB!  and  AgBs,  and  when  I = n/2  = cp.  From  a study  of  the  geometry 
considerations  (illustrated  in  Figure  10),  one  may  readily  determine  suitable 
sensor  locations  which  will  provide  tolerable  localization  resolution  over 
specified  areas  of  the  sphere. 
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For  convenience,  the  dependence  of  the  localization  resolution  area 
and  the  length  of  the  long  diagonal  on  the  angle  cp,  is  graphically  displayed 
i.i  Figure  11.  The  deterioration  of  localization  resolution,  for  cp  below 
about  ten  to  twenty  degrees,  is  quite  apparent  from  the  graph. 

6.0  EXAMPLES  OF  THE  USE  OF  THE  LOCALIZATION  ERROR  DATA 

To  understand  the  use  of  the  localization  error  data,  and  to  develop  a 
general  "feel"  for  the  magnitude  of  the  localization  resolution  which  is 
achievable  in  practice,  several  examples  will  be  given  in  the  areas  of 
electromagnetic  wave  propagation  and  in  underwater  acoustic  propagation. 

6.1  Electromagnetic  Localization  Applications 

In  the  following  applications,  it  will  be  assumed  that  electro- 
magnetic propagation,  which  essentially  follows  the  curvature  of  the  earth, 
is  employed  to  locate  a source  at  various  positions  with  respect  to 
strategically  located  receiving  stations.  The  two  pairs  of  receiving  sensors 
wi^l  be  assumed  to  each  consist  of  sensors  which  are  separated  on  great 
circle  paths  by  a distance  of  one-thousand  nautical  miles.  The  source  will 
be  located  at  various  distances  from  the  pairs  of  receiving  sensors  to  study 
the  resulting  localization  resolution  which  might  be  achieved  for  typical 
geometries . 

In  the  examples  to  be  presented,  we  shall  assume  a mean  global  radius 
of  3,440  nautical  miles  and  an  electromagnetic  propagation  velocity  around 
the  earth  surface  of  161,500  nautical  miles  per  second. 

Three  examples  of  source  ranges  (including  angles)  to  the  sensor  pair 
centers  (see  Figure  10)  are  postulated.  For  convenience,  all  of  the  pertinent 
data  on  the  source  resolution  are  tabulated  in  Table  2 for  easy  reference. 

The  upper  part  of  the  table  gives  the  postulated  data  an^  the  preliminary 
computations.  The  final  results,  illustrating  the  antit  ted  source 
resolution  parallelogram,  are  given  in  the  lower  five  rows  of  the  table. 

The  first  example  is  a near  field  problem  where  the  ranges  to  the  source 
are  less  than  the  separation  of  the  signal  sensors.  The  second  example  is  an 
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extremely  long  range  problem  where  the  source  is  located  in  the  region  of 
the  extreme  fa  field  convergence  zone.  In  the  third  example  the  source 
is  located  in  a region  where  the  localization  resolution  is  minimal.  The 
first  five  rows  of  the  table  postulate  the  general  geometry  of  the  problem. 

The  next  five  rows  ere  basic  data  which  are  computed  (or  derived  from  the 
earlier  graphs)  on  the  basis  of  the  individual  sensor  pair  configurations. 

The  final  five  lows  illustrate  the  pertinent  dimensions  of  the  resolutior 
error  parallelogram  which  are  computed  from  the  earlier  data.  (NOTE:  In 
the  subject  examples,  the  angle  cp  has  been  postulated  as  60°  in  all  cases. 

To  compute  this  angle,  additional  information  regarding  the  relative 
position  of  the  tuo  sensor  pairs  and  the  source  would  need  to  have  been 
specified. ) 

In  all  three  examples,  the  results  are  given  in  terms  of  ..he  time 
difference  resolution,  At,  (or,  alternatively,  2<tt)  measured  in  microseconds, 
for  a one  microsecond  time  difference  resolution,  then,  the  final  answer  is 
achieved  by  simply  letting  At  be  equal  to  unity.  For  other  time  difference 
resolutions,  the  final  answers  may  be  readily  derived,  it  is  significant  to 
note  that  the  localization  resolution,  for  examples  1 and  2 are  identical; 
in  spite  of  the  fact  that  the  relative  ranges  to  the  source  are  vastly 
different.  This  is  a result  of  the  convergent  zone  properties  of  the 
spherical  geometry.  Of  course,  the  subject  calculations  do  not  take  into 
account  propagation  anomalies  (which  must  appear  in  the  factor  At);  but 
even  so,  the  relative  accuracies  of  the  measurement  error  would  be  as  shown. 
This  means  that  although,  for  the  extremely  long  ranges  the  absolute 
accuracy  of  localization  may  be  somewhat  greater  than  shown  (due  to  uncertain- 
ties ir.  propagation  paths  and  velocities),  the  relative  position  of  two  (or 
more)  sources,  for  station  keeping  purposes,  would  be  as  normally  computed. 

6.2  Underwater  Acoustic  Localization  Applications 

In  addition  to  electromagnetic  propagation  applications,  the 
subject  theory  is  also  applicable  to  the  field  of  underwater  acoustics.  In 
these  latter  applications,  low  frequency  sound  will  travel  over  long  distances 
in  essentially  great  circle  paths  at  a propagation  velocity  of  approximately 
0.8  nautical  miles  per  second.  These  examples  of  underwater  acoustic  appli* 
cations  are  considered  in  the  same  manner  as  was  given  for  the  electromagnetic 
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:able  2 HESn.rS  OF  three  EM  localization  examples 


EXAMPLE 


#1 


1,000 
500 
75° 
400 
45  ° 


0. 114  AT 
37.5° 

0. 1 14  At 
18° 

0.013(At)‘ 


n 


1,000 

10,300 

75° 

10,400 

45° 


0, 114 AT 
37.5° 

0. 114 At 
18° 

0.013(At)? 


#3 


1 , 000, 

6,000 
20° 
4 , 000 
15° 


0.600, 

з. 5° 

0. 5 27  At 

и. 0° 

0.316<At)*| 


9 

o 

L; 

_3 

(h 

x. 


a 

Side  1 
Side  „ 
t(diag. ) 
Area 


(deg.) 
(NMi.) 
(NMi. ) 
(NMl. ) 
(Sq.  NMl.) 


60° 

0. 132A~ 

0.  3 32A'r 
0. 2 28 A7 
0.015(At): 


60° 

0.132AT 

0.132AT 

0.228AT 

().()15(At)? 


60° 
0.69 3 AT 
0.60'' At 
1.054A- 
0. 365(At)p 


NiTl’E-  A"  measure  ol  ambiguity  is  given  in  microseconds 


applications  discu..sed  in  the  previous  section.  The  pertinent  data  for 
these  examples  ate  displayed  in  Table  3. 

The  first  example  is  a near  field  problem  where  the  ranges  to  the 
signal  source  are  comparable  to  the  sensor  separation  distance,  0=200  NMi. 

In  tlte  second  example,  the  ranger  to  the  signal  source  have  bjen  increased 
to  ranges  of  approximately  1000  nautical  miles.  It  can,  therefore,  be 
anticipated  that  the  resulting  localization  tesolution  will  degrade 
appreciably.  In  the  third  example,  the  ranges  to  the  source  are  the  same 
as  for  example  2,  and  the  sensor  separation  distance  has  been  increased 
to  )U0  NMi.  This  last  example  Is  Intended  to  illustrate  the  improvement 
in  localization  resolution  which  can  be  realized  by  employing  more  widely 
separated  sensors  in  the  sensor  pair  coni  igur  at  ion In  all  of  the  underwater 
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TABLE  3 RESULTS  OF  THREE  UNDERWATER  ACOUSTIC  LOCALIZATION  EXAMPLES 


1 


EXAMPLE 

n 

n 

#3 

Di. 

d9 

(NMi. ) 

200 

200 

500 

Rl 

(NMi. ) 

300 

1200 

1200 

®i 

(deg.) 

20° 

15* 

15* 

R* 

(NMi.) 

250 

900 

900 

(deg.) 

30° 

20* 

20" 

*As 

i 

(NMi.) 

1 . 36At 

5.Q8AT 

2.04AT 

aui 

(deg.) 

18* 

14.5° 

14.5" 

As, 

(NMi.) 

1.20AT 

3.85AT 

1.54AT 

au? 

(deg.) 

26* 

19.5" 

19.5" 

A8j 

As, 

(Sr. NMi.) 

1.63(AT)a 

19.56(AT)S 

3. 14( At)3 

h 

CP 

(deg.) 

O' 

o 

o 

60* 

60" 

kJ 

Side  1 

(NMi.) 

1.57 At 

5.87AT 

2.36AT 

z 

o 

Side  2 

(NMi.) 

1.39At 

4.45AT 

1.78At 

§ 

t (di« t . ) 

(NMi.) 

2, 7 2 At 

8. 89 A. 

3.57At 

s 

JL 

Area 

(Sq.NMi.) 

1.88 (At) 3 

26. 6 (At) 3 

3.62(At)3 

* NOTE:'  measure  of  ambiguity  given  in  seconds 


acoustic  examples,  the  time  dlffeience  resolution  measure,  At,  is  given  in 
seconds. 

7 . 0 OPTIMUM  LOCALIZATION  RESOLUTION  SYSTEMS 

In  following  the  analysis  of  localization  resolution,  one  may  note  that 
there  exists  an  optimum  geometry  for  source  localization;  beyond  which,  the 
resolution  cannot  be  improved  by  geometric  considerations.  This  optimum 
geometry  is  realized  when  the  source  is  located  anywhere  on  the  great  circle 
path  Joining  the  two  sensor  locations.  On  this  line,  the  source  localization 
resolution  becomes,  simply, 


As  = (c/2)At 


(7-1) 
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As  the  source  deviates  from  this  line,  the  localization  resolution  will 
degrade  from  that  given  in  (7-1).  However,  the  amount  of  degradation  will 
be  dependent  on  the  great  circle  distance  between  the  two  sensor  locations 
and  the  relation  of  the  source  to  the  two  points.  It  can  reedil’-  be  shown 
th.it,  if  the  source  lies  within  a circle  which  passes  through  the  two  sensor 
points  and  is  orthogonal  to  tiie  great  circle  path  between  the  two  sensor 
points  (that  is,  the  straight  line  joining  the  two  sensor  points  forms  the 
diameter  of  the  circle),  then  the  localization  resolution  will  be  bounded 
by  the  following  relation, 

(c/.MA-  < As  < (c/./T)At  (7-2) 

ih i s tact  is  readily  evident  in  Figure  6.  Thus,  it  appears  that  for,  at 
least,  limited  regions  of  surveillance,  the  source  localization  resolution 
i an  he  optimized  ’ml  made  essentially  independent  of  source  range  and 
sensor  spearatlon  distance.  Two  applications  of  this  optimum  resolution 
phenomenon  will  now  be  considered. 

1 , 1 Channel  Navigation  Application 

Consider  a relatively  narrow  channel  through  which  vehicles;  such 
as,  ships,  aircraft,  etc.  may  wish  to  transit  under  controlled  navigation 
conditions.  By  placing  a pair  of  sensors  (or  sources)  across  and  orthogonal 
to  the  channel,  then,  the  optimum  in  localizacion  resolution  could  be 
achieved  foi  vehicles  passing  through  the  channel.  If  the  length  ot  the 
channel  Is  short  compared  to  the  channel  span,  one  pair  of  sensors  (sources) 
would  suffice  lor  the  navigation  system.  Otherwise,  appropriately  coded 
sets  of  sensors  could  be  placed  along  the  channel  to  guide  the  traffic  along 
the  entire  length.  The  lateral  localization  resolution  that  is  realizable 
by  such  a system  is  given  by  equation  (7-2).  For  example,  if  an  underwater 
acoustic  system  is  employed  across  a relatively  narrow  straight,  the  locali- 
zation resolution  within  the  straight  would  be, 

As  - 0,4A'r  (nautical  miles) 

800AT  (yards)  (l-i) 
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where  At  19  measured  in  seconds.  Thus,  if  the  measurement  of  time  register, 
t,  is  accurate  to  0.1  second,  the  cross-channel  resolution  would  be  80  yards. 
Shipping  lanes  (through  the  straight)  of  about  200  yards  in  width,  appear 
readily  assignable  for  the  control  of  local  traffic  in  this  application. 

Such  a system  might  be  feasible  tor  traffic  control  during  times  of  rela- 
tively poor  visibility  to  supplement  or  augment  electromagnetic  means  of 
navigation. 

7.2  Global  Navigation  Application 

A most  interesting  application  for  the  use  of  optimum  localization 
resolution  arises  in  the  case  where  the  sensor  (or  source)  locations  are 
spaced  one  hundred  and  eighty  degrees  apart,  on  opposite  sides  of  the  sphere. 
This  situation  was  described  in  Section  II,  3.1.1.  In  this  circumstance, 
the  source  (or  single  receiver)  will  always  lie  on  a great  circle  path  which 
joins  the  two  sensor  (or  source)  points.  Thus,  equation  (7-1)  represents 
the  localization  resolution  achievable  everywhere  on  the  «■  rface  of  the 
globe.  In  this  situation,  all  of  the  equi-time-register  curves  will  be 
circles  around  the  sphere,  similar  to  the  circles  of  constant  Latitude  on 
the  earth.  It  is  therefore  intriguing  to  consider  three  orthogonal  sets  of 
sensors  (sources),  located  on  orthogonal  axes  and  spaced  one  hundred  and 
eighty  degrees  about  the  sphere,  in  an  ideal  type  of  global  navigation 
scheme. 

As  an  example  of  the  proposed  global  navigation  system,  consider  that 
three  pairs  of  electromagnetic  transmitters  are  located  at  the  following 
earth  coordinates  relative  to  true  North. 

1.  (90°N,0°W)  and  (90®S,0°W)  True  North  and  True  South  Poles. 

2.  (0“N, 90°W)  and  (0°N,90oE)  Equator  at  90°West  and  90°  East  Longitude. 

3.  (0°N,0°W)  and  (0°N,180°W)  Equator  at  zero  and  180°  West  Longitude. 

The  frequency  of  the  transmissions  will  need  to  be  such  that  the  EM  propa- 
gation will  effectively  circle  the  earth  and  be  receivable  around  the  globe. 

For  the  hypothesized  system,  the  localization  equations  will  be,  simply, 
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60  !•  ujH/2  = sin'1!  J coss(uaTT/2;-sin2(u3Tr/2)  j 


= cos'1  j sin  (usn/2)+sin?(u5TT/2)'l 


(7-4) 


= sin'1 


= cos' 


1 


» tan 


where  ( 0O , 0P ) is  the  location  of  the  receiving  station  and, 

1 

l» i m cTj/np 
Up  - CTg/TTp 
U3  » C t3  /TTp 


(7-5) 


(7-5) 


The  subscripts,  in  the  above  Indicated  measures,  represent  the  three 
listed  transmitter  pairs  given  previously. 


For  the  postulated  navigation  system,  the  localization  resolution  in 
each  of  two  dimensions  will  be, 


As  * 0.081At 

- 162AT 


(nautical  miles) 
(yards) 


(7-7) 


where  L"  is  the  time  register  resolution  in  microseconds.  Thus,  for  about 
a ten  microsecond  time  register  resolution,  the  positional  resolution  of 
receiver  location  will  be  within  an  area  less  than  one  square  nautical  mile. 
Of  course,  due  to  KM  propagation  anomolies,  the  absolute  precision  of 
location  may  be  less  than  this.  However,  for  the  same  general  area,  the 
propagation  conditions  should  be  identical;  so  that,  the  above  localization 
resolution  Is  applicable  to  relative  position  control  for  station  keeping 
purposes,  etc. 
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V 


RESOLUTION  IN 
SOURCE  DYNAMICS 

In  Section  III,  we  developed  the  algorithms  or  relations  for  computing 
the  source  dynamics  (heading  and  speed)  from  a knowledge  of  the  system 
geometry  and  the  time  scale- factor  measures  between  two  pairs  of  receiving 
sensors.  Before  the  value  of  these  computations  can  be  assessed,  however, 
it  is  essential  to  know  the  confidence  or  trust  that  can  be  placed  in  the 
computed  results.  That  is,  how  accurate  are  the  resulting  dynamic  calcu- 
lations? If  we  can  assume  a negligible  geometry  and  computational  error, 
the  principal  contributor  to  error  (or  uncertainty  in  the  computed  results) 
will  be  the  accuracy  or  Inherent  resolution  of  the  basic  time  scale-factor 
measure,  6.  Thus,  what  we  intend  to  develop  in  this  section  is  the 
influence  of  errors  in  the  time  scale-factor  measure  on  the  computed  values 
for  the  source  speed  and  heading.  That  is,  we  shall  determine  the  resolu- 
tion in  the  measures  of  the  source  dynamics  in  terms  of  the  resolution  of 
the  basic  time  scale-factor  measures. 

The  problem  we  are  confronted  with  is  quite  complex  and  involved  since 
it  is  intuitively  evident  that  the  answers  we  are  seeking  will  be  dependent, 
not  only  on  the  rather  extensive  geometry  of  the  source  and  multisensor 
receiving  system,  but  also  on  the  interrelation  of  the  dynamic  variables 
themselves.  Never-the- less , by  a careful  selection  of  system  parameters 
which  convey  physical  meaning  within  the  geometric  system,  we  shall  be 
able  to  present  the  final  data  in  a form  readily  interpretable  and  useful 
in  physical  applications. 

1.0  BASIC  RELATIONS 

To  begin  our  solution,  we  shall  develop  the  basic  differential  equations 
which  relate  the  sensitivity  of  the  source  dynamics  to  variations  in  the 
time  scale-factor  measure. 

1. 1 Source  Velocity  Relation 

The  basic  differential  equation  for  the  source  velocity  will  be, 
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(1-1) 


dv 


5v 


5v 


cT6,d^  + av 


d 5, 


where,  in  the  present  stady,  the  subscript  1 will  be  employed  to  denote  the 
Lime  scale-factor  measure  for  the  first  pair  of  receiving  sensors,  and  the 
subscript  2 will  denote  the  time  scale-factor  measure  for  the  second  pair 
of  receiving  sensors.  Previously  (in  Section  III)  the  subscripts  for  these 
measures  were  xy  and  zw  respectively.  (For  convenience,  the  subscripts  1 
and  2 will  be  used  to  replace  the  earlier  subscript  notation  for  all  of 
the  system  parameters  dealing  with  the  two  pairs  of  receiving  sensors.) 


Using  equation  (3-7)  in  Section  III,  we  may  readily  carry  out  the 
differential  operations  called  for  in  (1-1)  to  obtain, 

dv  = c [b2  (ba5j  - bj  6aco82r|)d61  tb1  (b^  6a-b?  6,  cos2T))d6a] 

2b1bjsin2r|  J (b?  )2+(bj  ) ?-2bjba  5j  6acos2r) 


cP  fb2  (b7  -bj  6scos2'n)d61+b1  (bj  6g-bs  6j  cos2T))d  ) 
4 vbjba sinP2r| 


(1-2) 


1 • ^ Source  Heading  Relation 

Following  the  same  procedute  as  for  the  source  vel  'city,  the 
source  heading  differential  relation  may  be  derived  from  (3-6),  Section  III, 
as, 


d(3 


d6? 


-b1b?(6pd61-6,d6a)sin2n 

(b2  ) ?+(bj  6a )s-2bj b?  6?cos2r 


-cs‘(6gd61-61d68) 
4v2b;b?  sin2r, 


(1-3) 


2 . 0 K FLAT  IONS  FXPKFSSKD  IN  1NTKKPKKTABLK  FORM 

Although  tlie  basic  relations,  derived  above,  are  correct,  they  are  not 
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in  a form  that  is  readily  interpretable.  We  shall  therefore  proceed  to 
define  some  new  parameters  which  transform  the  parameters  given  in  (1-2)  and 
(1-3)  into  ones  which  are  readily  interpretable  in  terms  of  the  physical 
geometry  of  the  system  under  consideration. 


2.1  Definition  of  Geometric  Parameters 

In  Section  III,  we  defined  a number  of  parameters;  including, 

? » T1» ba,  which  were  convenient  in  developing  the  solution  for  the 
velocity,  v,  and  heading,  ft,  of  the  source.  We  shall  now  proceed  to  define 
a new  set  of  parameters  which  will  better  illustrate  the  geometry  of  the 
problem  we  are  dealing  with  and  provide  an  insight  into  the  effect  of  the 
system  geometry  on  the  resolution  of  the  measures  for  v and  9. 


Referring  to  equation  (3-3),  Section  III,  let, 

8l  * ? ‘ T|  • %(e0*-*-Boy  ) 
02  • ? + *1  “ fc(Bo*+Pow) 


(2-D 


The  above  two  new  parameters  represent  the  mean  angles  of  the  source 
rays  to  the  first  and  second  pairs  of  sensors  respectively. 


From  (2-1),  then,  it  is  seen  that, 

2r\  - 8s-8i  “ Aft?*  (2-2) 

where,  is  a new  parameter  representing  the  angle  between  the  two  mean 

angles  defined  in  (2-1), 


From  equation  (3-2),  Section  III,  we  may  write, 

bj  * s in% (60y 'Bo* ) * sin(A8i/2) 
ba  ■=  8in^(0ov-Bot ) = sin(Aea/2) 


(2-3) 
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where  Aflj  and  Ab9  represent,  respectively,  the  angle  between  the  rays 
emanating  from  the  source  to  the  two  sensors  in  each  pair  (that  is,  they 
represent  the  respective  pair  aperture  angles). 

It  will  also  be  convenient  to  define  one  additional  new  parameter  as 
follows.  Let, 


0ai  * + Bs>  * S (2-4) 

The  above  parameter  represents  the  mean  of  the  two  mean  angles 
Bj  and  g,  given  in  (2-1)  (which  is  also  the  mean  of  all  four  sensor  rays). 

To  obtain  a better  understanding  of  the  geometric  significance  of  the 
new  parameters,  they  are  presented  diagramatically  in  Figure  12.  (Although 
the  subject  diagram  is,  of  necessity,  drawn  on  a plane,  one  should  keep  in 
mind  that  the  lines  and  angles  are  intended  to  represent  what  would  be 
realized  geometrically  on  the  surface  of  a sphere.) 


For  a number  of  applications,  one  may  approximate  the  angles 
Agj/2  and  Ag3/2  by  the  following  equations, 

| Agj/2  ! tan-1  jD1coso,l /2R*  } 

! £B3/2  | tan-1  |Dacosa'.,/2R3  j 


(2-5) 


where,  Dj  and  D3  represent  the  great  circle  distance  between  the  two  sensors 
for  each  sensor  pair,  and  the  ranges  (Rj  and  R,)  and  angles  (»i  and  aa)  are 
as  shown  in  Figure  10.  The  suitability  of  the  above  approximations  will 
need  to  be  examined  for  the  specific  application  under  consideration. 


2. 2 Dynamic  Error  Relations 

We  may  now  employ  the  new  parameters,  defined  in  the  last  section, 
in  the  error  relations  given  in  (1-2)  and  (1-3).  To  do  this,  we  will  make 
use  of  equation  (3-4),  Section  III,  and  the  following  two  identities, 
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-sinA021co9(Bg-B) 


(2-6) 


sin(Bi -B)-sin(B?-B)cosA82i  = 

sin(B2-B)-sir.(e1-»)cosAesl  = sinAB21  cos  (b  5 - B ) * 

Carrying  out  the  resulting  algebra,  we  obtain  for  the  dynamic  error 
functions , 


(2-7) 


(2-8) 
(2-9) 

That  is,  the  variable  0 is  the  measure  of  the  source  heading  relative  to 
the  mean  angle  of  all  rays  from  the  source  to  the  four  receiving  sensors 
(see  Figure  12). 

(NOTE;.  The  heading  error,  dB , given  in  ( 2- 8 > is  in  radians.  To  convert  to 
degrees,  this  relation  should  be  multiplied  by  the  factor  180/n. ) 

2.3  Properties  cf  the  Dynamic  Error 

From  (2-7)  and  (2-8)  it  is  evident  that, 

dg  = -dv/v  (2-10) 

for  the  heading  angle  advanced  ninety  degrees.  It  is  further  evident  that 
the  resolution  in  source  heading  will  be  inversely  proportional  to  the 
source  velocity.  This  is  to  be  expected  intuitively;  since,  for  a source 
velocity  of  zero,  the  heading  becomes  completely  ambiguous. 


dv  = 


cost 


>Cfla±)  ■ d6,  . Sgg.Xft.i-il.  d6 

2sinAjjsl  L sin(ABi/2)  sin(ABg/2)  s 


and , 


f cos(0-ABqt  /2)  _ cos,(>AB~aI/2)  1 

L sin(Afl , / 2)  61  sin( ABj/2)  7 J 


2sinAgai  L sin(AB1/2) 


de  „ £-.—  [ sln,(L-Cg-)  . dfi  . sinXiirll..  d. 

2vsinA??i L sin(AB1/2)  • sin(ABs/2)  a 


where, 


r gjiiX3-^.i7.?).  d6  . d5 

2v8inAf?1l  sin(ABi/2)  1 sin(ABs/2)  a 


n = B-8al 
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It  Is  further  evident  that  for  some  particular  source  heading,  the 
dynamic  error  will  be  inversely  proportional  to  the  sine  of  the  angle  A@sl. 
Consequently,  the  dynamic  error  tends  to  minimize  when  gj  is  orthogonal 
to  02  (see  Figure  12).  (NOTE:,  The  fact  that  this  is  not  true  for  all 
source  headings  will  be  discussed  later.  At  this  point  we  are  confining 
our  remarks  only  to  the  general  trend  where  all  heading  angles  need  to  be 
considered.)  The  error  also  minimizes  as  sin^!/'?,)  and  sin(AB?/2) 
approaches  unity.  Thus,  the  error  reduces  as  the  sensor  separation  distance 
Dj  D?  increases  and  the  range  to  the  source  (from  the  sensor  midpoints) 
uecreases  (see  Figures  10  and  12).  These  results  appear  Intuitively 
acceptable  based  "pon  geometric  considerations. 

3.0  DYNAMIC  ERROR  MEASURE  STATISTICS 

Equations  (2-7)  and  (2-8)  can  be  used  to  determine  the  error  in  the 
source  dynamics  for  specified  error  in  the  measures  of  ft5  and  6S  (providing 
the  errors  are  quite  small.)  Of  greater  significance  will  be  the  expected 
values  of  these  errors  and  their  standard  deviations.  These  statistical 
measures  of  error  will  now  be  computed  in  terms  of  the  standard  deviation 
of  the  measure  of  the  time  scale-factor  shift  variable  ft. 

3.1  Probability  Density  Function  for  the  Dynamic  Error  Measures 

From  (2-7)  and  (2-8)  it  is  evident  that  the  dynamic  error  measures, 
dv  and  dg , take  the  form. 


f = A Aft}  + BAftj  (3-1) 

where  f represents  the  dynamic  error  statistic  (dv  or  dB),  A and  B are 
parameters,  and  Aftj  and  A6S  are  independent  statistical  samples  of  the  time 
scale-factor  shift  error  statistic.  We  wish  to  determine  the  mean  or 
expected  value  of  f and  also  the  variance  of  f about  this  expected  value. 

To  accomplish  our  objective,  we  shall  assume  that  the  time  scale-factor 
error  statistics,  A6j  and  A5S , are  independent  and  have  a zero  mean  Gaussian 
distribution.  For  this  assumption  (which  is  quite  reasonable),  it  is  a 
simple  matter  to  show  that  the  dynamic  error  statistic,  €,  also  has  a zero 
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mean  Gaussian  distribution.  The  expected  value  of  C-  will,  therefore,  be 
zero,  and  it  may  readily  be  shown  that  the  variance  of  C will  be, 


= (ae+bs)o63 


where,  H , is  the  standard  deviation  of  the  time  scale-factor  shift 
variable 


3.2  Standard  Deviation  of  the  Dynamic  Error  Measures 

Using  the  results  of  (3-2)  in  (2-7)  and  (2-8)  we  obtain  the 
standard  deviation  of  the  dynamic  errors  as, 


I cos2r-A0al/2)  cos2(0f  Aflsi  /2) 


2 I sinAg?1  | y sina(A8j/2) 


sins(AB?/2) 


'sin2(r-A8si/2)  sina(OfAial/2) 
+ 


nv  1 sinAe?x  ! \/  sinB(Ae1/2)  sin3(A8s/2) 


where  the  standard  deviation  of  the  source  heading  is  given  in  degrees. 

3 . 3 Properties  of  the  Standard  Deviations 

The  following  properties  of  the  standard  deviations  given  in 
(3-3)  and  (3-4)  will  be  found  useful. 


CJj  ('1;  Ajj?i  ,Agi , Aea) 


r-90;AB-sl,ABl,A8.,) 


'’v  (*  » Ap  a j , A8  j , A8  a ) — ( 1 8 0 - . , j A8  a ^ , A8  a ) 


( i-ASajjAB^jASg)  — ^ i Ab  a x , AB  a , AB  x ) 
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°v  (^;  180- , Aflj , Ab3  ) = f7v  (0-90;  Ags  j , Agg , Agj ) 


- oy  (90-0;  Aggi » Ag i > Ab3 ) 


(3-7) 


Equation  (3-5)  Illustrates  the  simple  relationship  between  the 
standard  deviations  of  the  two  dynamic  measures,  and  equations  (3-6)  and 
(3-7)  demonstrate  the  symmetries  that  exist  between  the  modified  heading 
variable,  0,  and  the  parameters  A83i,A8i  and  Ag3. 

3.3.1  Non-Dimensional  Equation  Form;  For  convenience,  the  following 
nondimen sional  form  will  be  used  to  express  the  source  dynamics  errors. 
Let, 


Cy  (^;  Agal  ,b) 


8in(Ae« /2) 


£v 

CCTfi 


and, 


4^/cos3(0/-A8ai/2)+b2cos3(0,+AB31/2)' 
2 ! sin(ABai) ! 


vnB 

sin(AC./2)  ! — 


180 

TT 


%(0-90;Ae?1,b) 


(3-8) 


(3-9) 


where. 


= ^lnLAS^j.2) 
sin(ASb/2) 


! A?.  | * | ABb 

and  where, 


(3-10) 


As,  is  the  smaller  of  the  absolute  values  of  Aflj  and  A03 
Asb  is  the  larger  of  the  absolute  values  of  Ab^  and  A83 
is  a new  modified  source  heading  angle  related  to  the 
modified  source  heading  angle  Q in  a manner  defined  in 
Table  5. 
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The  above  functions  are  plotted  graphically  in  Figures  13  through  17 
for  convenience  in  data  abstraction.  Selection  of  the  appropriate  scales 
(on  the  graphs),  as  a function  of  the  variables  and  the  range  on  the 
pertinent  parameters,  is  given  in  Table  4 for  ready  reference. 


FAB  Lb  4 INTERPRETATION  OF  SCALES  ON  THE  ACCOMPANYING  GRAPHS 


\ 

%*  ]sin(Ae, /2)  | CTv/(c06) 
(Use  left  hand  ordinate  scale) 

rB  = |sin(Ag. /2)|  vCs/(cOg) 
(Use  right  hand  ordinate 
scale) 

|Ab  i ! s 1 ABj  ' 

|abj  ! < ! Ae i! 

1 5 |ab?| 

!ab?  ! <-  |Aea  ! 

(.'e.-ABi) 

(Ae,=ABP) 

(Ab»=Ab1  ) 

(Ab.  = AB3) 

o 

VI 

... 

o'  ■ 1. 

(lower  scale) 

r\  f _ _ ***> 

(lower  scale) 

o'. 90  „ - 
(upper  scale) 

o'-90=  -o 
(upper  scale) 

H 

fV 

_ .A 

q'  » O 

f)'-90  - -C, 

o'. 90  E o 

1 <U 
•1 

Afig l < 0 

(lower  scale) 

(lower  scale) 

(upper  scale) 

(upper  scale) 

r '-90= 

o'-90  » 0 

o _o 

n ' _ o 

1 Ui 

(upper  scale) 

(upper  scale) 

(lower  scale) 

(lower  scale) 

VI 

r>  /-90  * 0 

o'- 90  = -D 

o ' _ 0 

- 

O 

CT> 

A6,i<  0 

(upper  scale) 

(upper  scale) 

(lower  scale) 

(lower  scale) 

NOTE:  See  Figure  12  for  geometric  interpretation  of  variable  and  parameters 


J . 1.2  Analysis  of  the  Graphical  Data:'  A study  of  the  standard 
deviations  displayed  graphically  in  Figures  13  through  17  reveals  the 
following:; 

(1)  The  dynamic  error  (or  standard  deviation)  will  be  directly 
proportional  to  the  sine  of  one-half  the  aperture  angle  of 
the  smallest  sensor-pair  aperture  (Ap»).  (The  aperture 
angle  of  a sens  r pair  is  taken  as  the  angle  between  the 
two  great  circle  rays  emanating  from  the  source  to  the  two 
individual  sensors  forming  the  sensor  pair.)  Therefore, 
the  larger  that  the  smallest  of  the  two  aperture  angles 
becomes,  the  better  will  be  the  resolution  of  the  dynamic 
measures.  The  optimum  aperture  angle  will,  of  course,  be 
180  degrees.  This  aperture  angle  occurs  when  the  source 
is  located  on  the  great  circle  path  between  the  sensor  pair. 
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(2)  The  dynamic  error  will  be  relatively  insensitive  to  the 
parameter  b (the  ratio  of  the  sines  of  the  two  sensor 
pair  aperture  half  angles).  Given  the  fact  that  the 
minimum  of  the  two  aperture  angles  is  fixed,  the  fact 
that  the  larger  of  the  two  aperture  angles  is  about 
equal  to  the  smaller,  or  is  very  much  larger,  is  of 
only  minor  significance.  Naturally,  the  greatei  the 
difference,  the  better  will  be  the  source  dynamics 
resolution.  But,  the  improvement  will  be  relatively 
small. 


(3)  The  dynamic  error  sensitivity  to  the  dual  pair  aperture 
angle,  A§P1,  will  be  strongly  dependent  on  the  source 
heading  angle  relative  to  the  mean  sensor  angle  gal 
(see  Figure  12).  Over  most  heading  angles,  the  dynamic 
error  increases  as  the  dual  pair  aperture  angle  becomes 
smaller,  but  the  increase  in  error  does  not  become 
excessive  until  this  dual  pair  aperture  angle  reduces 
to  something  less  than  about  fifteen  degrees.  For  dual 
pair  aperture  angles  below  about  five  or  six  degrees, 
the  loss  in  source  dynamics  resolution  can  become 
ohibitivnly  large.  An  interesting  and  unexpected 
result  is  that,  over  a particularly  narrow  range  of 
heaci,.lt,  angles,  the  source  dynamic  resolution  will  be 
relatively  independent  of  the  dual  pair  aperture  angle 
A8?1,  and  can  even  improve  with  a decrease  in  Agal. 

The  overall  effect,  though,  is  that  a decrease  in  the 
dual  pair  aperture  angle  will  result  in  a decrease  in 
the  source  dynamic  resolution  (increased  dynamic 
measure  error). 


(4)  It  was  mentioned  above  that,  for  dual  pair  aperture  angles, 
different  from  ninety  degrees,  the  resolution  of  the 
dynamic  measures  will  be  dependent  on  the  source  heading. 

An  interesting  facet  of  this  phenomenon  is  that  the 
velocity  error  and  the  heading  error  are  in  quadrature  in 
this  regard.  Hiat  is,  for  those  headings  which  minimize 
the  velocity  error,  the  heading  e.rror  will  be  maximized. 
And,  conversely,  where  the  heading  error  is  a minimum,  the 
velocity  error  will  be  maximum.  Thus,  there  appears  to  be 
a trade-off  in  resolution  between  the  vt locity  measure  and 
the  heading  (.or  course)  measure  as  a function  of  the  source 
heading.  Naturally,  there  exists  source  headings  which 
result  in  about  equal  resolution  of  both  dynamic  measures 
(where  neither  is  maximum  or  minimum). 

Cj)  Another  interesting  phenomenon  is  that  when  the  absolute 
value  of  the  dual  pair  aperture  angle,  )g?1,  is  less  than 
ninety  degrees,  the  minimum  error  for  source  velocity  is 
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realized  for  source  headings  which  are  normal  to  the 
mean  sensor  direction,  g?1  (see  Figure  12).  Thus,  the 
minimum  error  for  source  heading  will  occur  for  source 
headings  along  the  mean  sensor  direction,  gal.  On  the 
other  hand,  when  the  absolute  value  of  the  dual  pair 
aperture  angle,  Agal,  becomes  greater  than  ninety 
degrees,  the  converse  will  be  true.  That  is,  in  this 
circumstance,  the  minimum  error  for  source  velocity 
will  occur  for  source  headings  along  the  mean  sensor 
direction,  5ai,  and  the  minimum  error  for  source 
heading  will  occur  for  source  headings  orthogonal  to 
this  direction. 


(6)  As  an  overall  observation,  it  can  be  concluded  that  the 
source  dynamics  resolution  will  be  near  optimum  when  the 
two  sensor  pair  aperture  angles,  Ab3  and  Aga,  are  large 
(preferably  near  180  degrees)  and  the  dual  pair  aperture 
angle,  A0al,  is  in  the  neighborhood  of  90  degrees.  As  a 
consequence,  when  the  source  is  located  in  the  near  field 
of  the  sensors;  so  as,  to  be  effectively  surrounded  by  the 
sensors,  the  resolution  of  the  dynamic  measures  will  be 
high.  As  the  source  location  becomes  remote  from  the 
general  sensor  configuration,  Lhe  source  dynamic  resolution 
can  be  expected  to  decrease.  But,  this  decrease  in  source 
dynamic  resolution  will  continue  only  to  the  point  where 
the  source  is  removed  approximately  one-quarter  of  the  way 
around  the  spherical  globe  from  the  geometric  center  of  the 
sensor  configuration.  As  the  source  location  becomes 
further  removed,  the  source  dynamic  resolution  can  be 
expected  to  improve  again.  And,  as  the  source  reaches  the 
extreme  far  field,  on  the  opposite  side  of  the  globe  from 
the  sensor  configuration,  the  relevant  aperture  angles  will, 
once  again,  achieve  rhe  same  values  they  took  on  in  the  near 
field.  Thus,  the  convergent  zone  properties  discussed  in 
regard  to  the  sensor  location  resolution,  will  also  apply 
to  the  sensor  dynamics  resolution.  Consequently,  it  appears 
possible  to  achieve  an  unbelievably  high  heading  and  speed 
resolution  for  a source  located  at  extremely  long  great 
circle  ranges  from  the  relevant  sensor  configuration  when 
we  are  working  in  a spherical  geometry.  The  practicality  of 
this  seeming  paradox  will  be  dependent  on  the  ability  to 
realize  reliable  signal  propagation  over  the  extremely 
long  propagation  paths. 


4 . 0 EXAMPLES  OF  THF.  USK  OF  THE  DYNAMICS  KKHOR  DATA 

To  understand  the  use  of  the  source  dynamics  error  data,  and  to  develop 
a general  "feel"  for  the  ignitude  of  the  velocity  and  heading  resolution 
which  is  achievable  in  practice,  several  examples  will  be  given  in  the  areas 
of  electromagnetic  wave  propagation  and  in  underwater  acoustic  propagation. 
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4.1  electromagnetic  Source  Dynamics  Application 

In  the  following  applications,  it  will  be  assumed  that  electro- 
magnetic propagation,  which  essentially  follows  the  earth's  curvature,  is 
employed  to  measure  the  source  dynamics  (heading  and  speed)  of  an  aircraft 
in  motion.  The  location  and  dynamics  of  the  aircraft  will  be  varied  with 
respect  to  strategically  located  receiving  stations.  The  geometry  of  the 
receiving  stations  will  also  be  varied  to  demonstrate  the  effect  of  these 
parameters . 

Four  examples  have  been  chosen  for  study,  and  the  results  are  tabulated 
in  Table  5 for  easy  reference.  The  upper  portion  of  the  table  gives  the 
postulated  geometry  and  source  dynamics.  The  central  position  of  the  table 
lists  the  results  of  the  preliminary  computations.  And,  the  lower  portion 
of  the  table  illustrates  the  final  results  for  the  resolution  of  the  target 
dynamics.  The  ! Mndard  deviation  of  the  target  dynamics  is  first  given  in 
terms  of  the  standard  deviation  for  the  basic  time  scale-factor  measure, 
and  then  followed  by  actual  numbers  based  upon  an  assumed  value  for  ccr^ 
which  appears  reasonable  to  achieve. 

The  first  example  (see  Table  5)  depicts  a near  field  situation  where 
the  ranger  to  the  source  are  less  than  the  separation  distance  (Dj  and  Ds) 
between  sensor  stations.  The  source  velocity  (speed  of  the  aircraft)  is 
given  as  500  KTS.  The  second  problem  is  identical  to  the  first,  except 
that  the  source  has  been  translated  to  the  extreme  far  convergence  zone 
field  on  the  opposite  side  of  the  globe.  It  will  he  noted  that  the  resolu- 
tion of  the  source  (aircraft)  dynamics  is  the  same  in  both  cases. 

In  the  third  and  fourth  examples,  the  sensor  separation  distance  has 
been  reduced  to  two  hundred  nautical  miles.  In  the  third  example,  the 
source  ranges  can  be  considered  as  intermediate,  and  in  the  fourth  example 
they  would  be  considered  as  reasonably  long.  The  standard  deviation  of  the 
source  dynamics  can  he  expected  to  degrade  measurably  due  to  the  decrease 
in  the  source  pair  aperture  angles.  Th's  is  moderated,  somewhat,  by  the 
fact  that  the  dual  pair  aperture  angle  has  been  increased  to  ninety  degrees. 
In  addition,  the  higher  source  speed  (600  KTS  as  compared  with  250  KTS)  nas 
resulted  in  a higher  resolution  for  the  source  heading. 
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TABLE  5 RESULTS  OF  FOUR  EM  DYNAMIC  EXAMPLES 


EXAMPLE 

n 

#2 

#3 

#4 

Dj  ,Dj 

(NMi) 

1,000 

1,000 

200 

200 

R1 

(NMi) 

500 

10,300 

400 

1,200 

(deg.) 

30° 

o 

O 

CO 

18° 

5° 

(NMi) 

400 

10,400 

300 

800 

(deg.) 

45° 

45° 

30° 

15° 

V 

(KTS.) 

500 

500 

250 

600 

(deg.) 

40° 

40° 

60° 

20° 

ABj 

(deg.) 

82° 

82° 

27° 

mm 

ASs 

(deg.) 

83° 

83° 

32° 

13.8° 

sln( A8x/2) 

0.66 

0.66 

0.233 

0.083 

sin(A6a/2) 

0.66 

0.66 

0.277 

0.120 

b? 

1 

1 

-- J 

0.707 

0.478 

r* 

tv 

3 

(deg.) 

45° 

45° 

30° 

90° 

Oy 

(KTS) 

1.14(coft) 

1.14(cct?,) 

3.22(Cog) 

5.78(caj) 

(deg. ) 

0. 134(co*) 

0.114(cO6) 

1.06(coj) 

0.437(co6) 

(c  C6) 

?v  (KTS) 

11.4  KTS 

11.4  KTS 

32.2  KTS 

57.8  KTS 

= 10  KTS 

(deg.) 

1.14° 

1.14° 

o 

O' 

o 

4.37" 

The  above  examples  illustrate  that  quite  useful  results  should  be 
achievable  over  a variety  of  system  geometries  and  source  dynamics. 

4.2  Underwater  Acoustic  Applications 

Since  the  subject  theory  is  applicable  to  underwater  acoustic 
propagation,  we  have  considered  four  examples  in  this  area  for  demonstration. 
The  examples  are  presented  in  tin  same  manner  as  for  the  electromagnetic 
propagation  applications  described  in  the  last  section.  The  resulting  data 
are  tabulated  in  Table  6 for  convenient  reference.  (In  the  case  of  the 
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TABLE  6 RESULTS  OF  FOUR  UNDERWATER  ACOUSTIC  DYNAMIC  EXAMPLES 


EXAMPLES 


#1 

#2 

#3 

#4 

Da,D?  (NMi, ) 

180 

180 

200 

500 

R,  (NMi.) 

200 

380 

1200 

1200 

(deg.) 

30° 

13° 

15° 

15“ 

(NMi.) 

250 

400 

900 

900 

n?  (deg.) 

20° 

10“ 

20“ 

20“ 

v (KTS. ) 

10 

20 

25 

25 

' (deg.) 

' lJ  ' ‘ 
20° 

45° 

-30° 

- 

-30“ 

4i  (deg.) 

42° 

26c 

9.2“ 

22.8° 

37° 

25° 

12.0“ 

29.2° 

s in  ( Ae  i / 2 ) 

0.363 

0.225 

0.080 

0.198 

- - 

sin(AB->/2) 

0.321 

0.216 

0.105 

.. 

0.252 

b' 

0.78 

0.925 

0.589 

0.617 

4a  i (deg.) 

40° 

53° 

60“ 

60° 

Uv  (KTS.) 

2.96(cr-f;) 

2.92(ca5) 

5 .63  (ctr«s) 

2.27(cr6) 

4 (deg.) 

11.6(ca6) 

11. 2(cn^) 

14 . 9 (cCT5) 

6.03(ca6) 

(c^r) 

= 0.15  Kl’S. 

f'v  (KTS.) 

0.44  KTS. 

0.44  KTS. 

0.84  KTS. 

0.34  KTS. 

(deg.) 

1.74° 

1.68° 

2.24° 

0.90“ 

underwater  acoustic  applications,  a value  for  c^ft  of  0,15  was  assumed  to  be 
realizable  and  is  employed  in  the  last  two  rows  of  the  table.) 

'n  the  first  two  examples  a sensor  separation  distance  (D1  and  D?)  of  180 
nautical  miles  was  cho.,«n  and  source  ranges  in  the  near  to  intermediate  tield 
are  presumed  l’he  resulting  ->t  mdard  deviation  for  the  dynamic  variables 

84 


(source  speed  and  heading'*  are  comparable  and  fall  within  a range  which 
appears  quite  attractive. 

In  the  latter  two  examples,  the  source  range  has  been  increased  to  the 
neighborhood  of  one  thousand  nautical  miles.  The  two  examples  are  identi- 
cal with  the  exception  that  the  sensor  separation  distances  'Dj  and  Da)  have 
been  increased  from  200  nautical  miles  (in  example  #3)  to  500  nautical  miles 
(in  example  #4).  The  results  demonstrate  the  improvement  achieved  by  the 
wider  sensor  separation  at  the  longer  ranges. 

The  above  examples  demonstrate  the  usefulness  of  time  scale- factor 
measure  to  determine  source  dynamics  in  underwater  acoustic  applications, 
providing  the  resolution  of  the  time  scale-factor  measure,  6,  can  be  made 
sufficiently  high. 

5.0  OPTIMUM  SOURCE  DYNAMICS  RESOLUTION  SYSTEMS 

In  following  the  analysis  of  source  dynamic  resolution,  one  may  note 
that  there  exists  a near  optimum  geometry  for  source  dynamics  resolution; 
beyond  which,  the  resolution  cannot  be  improved  by  geometric  considerations. 
This  optimum  geometry  is  realized  when  Ag1  and  Aga  are  both  180  degrees, 
and  Is  90  degrees  (see  Figure  12).  When  this  situation  exists, 

equations  (3-3)  and  (3-4)  reduce  simply  to, 

<?v  = ^(cn6)  (5-1) 

and , 

90/  _ . 28.65,  _ , ,r 

°8  “ ^<c r,6)  = — — ($~2) 

The  above  relations  will  be  realized  when  the  source  is  located  on 
the  great  circle  paths  which  Join  the  sensors  in  both  pairs  of  receiving 
stations  and  when  these  great  circle  paths  are  orthogonal.  When  the 
latter  requirement  is  not  fully  met,  the  result  will  not  be  seriously 
changed  providing  A6?i  does  not  become  significantly  less  than  about  45 
degrees  (see  Figure  13).  As  a consequence,  relations  (5-1)  and  (5-2)  will 
be  applicable  for  the  global  navigation  system  described  in  Section  IV,  7.2. 
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For  a more  localized  sicuation,  where  a particular  area  is  to  be 
covered,  the  optimum  sensor  geometry  would  be  as  illustrated  in 
Figure  18. 
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Fig.  18  Optimum  Sensor  Geometry  for  a Localized  Area  Coverage 

In  Figure  18  the  sensor  pairs  are  xy  and  zw,  and  the  desired  area  of 
coverage  lies  within  the  circle  joining  all  of  the  sensors.  Naturally,  the 
closer  the  source  is  to  the  center  of  the  circle,  the  more  nearly  optimum 
will  be  the  resolution.  However,  it  can  be  shown  that  for  the  source  located 
anywhere  within  the  indicated  circle,  the  standard  deviation  of  the  velocity 
error  will  be  bounded  by  the  relation, 

0.5(cC*)  ^ oY  < 1.3(cCTj)  (5-3) 

As  it  turns  out,  the  same  general  system  geometry  which  optimizes  the 
source  localization  resolution  will  also  be  optimum  for  the  source  dynamics 
resolution. 
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